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PKEFACE. 



The object of the following treatise is to exhibit, in a cone 
form, the elementary properties of the expressions known 
the name of Laplace's Functions, or Spherical Harmoni 
and of some other expressions of a similar nature. I do n 
of course, profess to have produced a complete treatise 
these functions, but merely to have given such an introdi 
tory sketch as may facilitate the study of the numerc 
works and memoirs in which they are employed. 
Spherical Harmonics derive their chief interest and util 
from their physical applications, I have endeavoured fn 
the outset to keep these applications in view, 

I must express my acknowledgments to the Eev. C. 
Prior, Fellow of Pembroke College, for his kind revision 
the proof-sheets as they passed through the press. 



N. M. FERRERS, 



QovnujR 4ND Oaius Gollbgb, 
August, 1877. 



F. H. 
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CHAPTER I. 

INTRODUCTORY. DEFINITION OP SPHERICAL HARMONICS. 

1. If V be the potential of an attracting mass, at any 
point Xy y, z, not forming a part of the mass itself, it is 
known that Fmust satisfy the differential equation 

daf^^^ d^ ^ ^^J' 

or, as we shall write it for shortness, V*F=0. 

The general solution of this equation cannot be obtained 
in finite terms. We can, however, determine an expression 
which we shall call V^, an homogeneous function of x, y, z 
of the degree t, i being any positive integer, which will 
satisfy the equation ; and we may prove that to every such 
solution FJ there corresponds another, of the degree — (i + 1), 

expressed by -j^ , where r* = a?* + y' + z\ 

For the equation (1) when transformed to polar co-ordi- 
nates by writing a; = r sin ^ cos ^, y = r sin ^ sin (p, z^r cos 0, 
becomes 

And since V satisfies this equation, and is an homo- 
geneous function of the degree i, F, must satisfy the equa- 
tioi^ 

... ,, ^. 1 d f\ ^dVA . 1 d'V, ^ 
^^*HmddO\ d$ J ^md a^'' 



F. H. 
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INTKODUCTORT. 



since this is the form which equation (2) assumes when V^ 
is an homogeneous function of the degree u 

Now, put F, = r^*** PJ, and this becomes 
.•(.■+l)H-.R+3l-,|(H-.h.»f)+jij|.(^R).«. 



or 



Now, since Z7] is a homogeneous function of the degree 

-(»■ + !). 

therefore equation (2) becomes 



ar" 



^ 1 rf /. ^dU,\^ 1 ^t^_o 

■^sli^^ V^'''^:^;7"sin''e^ d<i? "■^' 



shewing that 17, is an admissible value of F, as satisfying 
equation (2)^ 

It appears therefore that every form of U^ can be ob- 
tained from F„ by dividing by r'**\ and conversely, that 
every form of V^ can be obtained from U^ by multiplying 
by r*""\ Such an expression as F] we shall call a Solid 
Spherical Harmonic of the degree i. The result obtained 
by dividing F, by t\ which will be a function of two inde- 
pendent variables 6 and ^ only, we shall call a Surface 
Spherical Harmonic of the same degree. A very important 
form of spherical harmonics is that which is independent 
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DEFINITION OF SPHERICAL HARMONICS, S 

of ^. The solid harmonics of this form will involve two of 
the variables, x and y, only in the form oj* + ^, or will be 
functions of a?* + y* and z. Harmonics independent of ^ are 
called Zonal Harmonics, and are distinguished, like spherical 
harmonics in general, into Solid and Surface Harmonics. 
The investigation of their properties will be the subject of 
the following chapter. 

The name of Spherical Harmonics was first applied to 
these functions by Sir W. Thomson and Professor Tait, in 
their Treatise on Natural Philosophy. The name " Laplace's 
Coefficients" was employed by Whewell, on account of 
Laplace having discussed their properties, and employed 
them largely in the Mecanique Vdeste. Pratt, in his 
Treatise on the Figure of the Earth, limits the name of 
Laplace's Coefficients to Zonal Harmonics, and designates 
all other spherical harmonics by the name of Laplace's 
Functions, The Zonal Harmonic in the case which we shall 
consider in the following chapter, i.e., in which t3ie system 
is symmetrical about the line from which is measured, 
was really, however, first introduced by Legendre, although 
the properties of spherical harmonics in gieneral were first 
discussed by Laplace; and Mr Todhunter, in his Treatise, 
on this account calls them by the name of "Legendre's 
Coefficients," applying the name of "Laplace's Coefficients" 
to the form which the Zonal Harmonic assumes when in 
place of cos 6, we write cos 6 cos ff + sin 6 sin 6' cos (^ — ^'). 
The name " Kugelfunctionen " is employed by Heine, 
in his standard treatise on these functions, to designate 
Spherical Harmonics in general. 



1—2 
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CHAPTER 11. 

ZONAL HARMONICS. 

1. We sliall in this chapter regard a Zonal Solid Har- 
monic, of the degree i, as a homogeneous function of 
(«* +y^ii and z, of the degree i, which satisfies the equation 

d'V , d^V d'r_^ 

da?^ df^H^"^' 

Now, if this be transformed to polar co-ordinates, by 
writing r sin 6 cos (f> for a?, r sin sin ^ for y, rcosO for z, we 
observe, in the first place, that a;* + y' = r* sin' ^. Hence 
V will be independent of <f}, or will be a function of r 
and only. The differential equation between r and 
which it must therefore satisfy will be 

dr sm^ dd \ d0j 

Now V, being a function of r of the degree i, may be 
expressed in the form r*F^, where P, is a function of only. 
Hence this equation becomes 

or, putting cos = fi, 

In accordance with our definition of spherical surface 
harmonics, P, will be the zonal surface harmonic of the 
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ZONAL HAKMONICS. 5 

degree t. When it is necessary to particularise the variable 
involved in it, we shall write it P, iji)^ 

The line from which is measured, or in other words 
for which /* = 1, is called the Axis of the system of Zonal 
Harmonics; and the point in which the positive direction 
of the axis meets a sphere whose centre is the origin of 
co-ordinates, and radius unity, is called the Pole of the 
system. 

Any constant multiple of a zonal harmonic (sojid or 
surface) is itself a zonal harmonic of the same order. 

2, The zonal harmonic of the degree i, of which the 
line fi = 1 is the axis, is a perfectly determinate function of 
fi, having nothing arbitrary but this constant. For the 
expression r*P^ may be expressed as a rational integral 
homogeneous function, of r and z, and therefore P, will be 
a rational integral function of cos 0, that is of /a, of the 
degree i, and will involve none but positive integral powers 

of /Lt. 

But P, is a particular integral of the equation 

and the most general form of f(fi) must involve two ar- 
bitrary constants. Suppose then that the most general 

form of /(/a) is represented by P^ jvdfjL. We then have 

+ 2(l-;.«)g. + P,|{(l-;.')..}. 

Hence, adding these two equations together, and ob- 
serving that, sincQ P, satisfies the equation (3), the coefficient 
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fr ZONAL HARMONICS. 

of jvdfM will be identically equal to 0, we obtain, for the de* 
tiennination of v, the equation 

r.|{(i-/*')«}-h2(i-;.')^'.=o. 

whence -P,(1-m')^+2|(1-/**)^-/*p]» = 0, 

the integral of which is 

log V + log P* (1 — /**) = log Cj = a constant ; 
0. 



t>= ^^L 



Hence Jed/* = (7+ O; j p^. ^^^^i^ J 

and we obtain, for the most general form oif(ji), 

Now, P| being a rational integral function of /i of i 

dimensions, it may be seen that /t^- — rrpi ^^^ assume the 

form of the sum of i + 2 logarithms and i fractions, and 
therefore cannot be expressed as a rational integral function 

of /i. Expressions of the foi:m P, / . ^ pa ^'^^ called KugeU 

functionen der zweiter Art by Heine, who has investigated 
their properties at great length. They have, as will hereafter 
be seen, interesting applications to the attraction of a sphe- 
roid on an external point. We shall discuss their properties 
more fully hereafter. 

3, We have thus shewn that the most general solution 
of equation (2) of the form of a rational integral fi^nction of u. 
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toNAL fiARJUOis^ICS; 7 

involves but one arbitrary constant, and that as a factor. 
We shall henceforth denote by P,, or P, (ft), that particular 
form of the integral which assumes the value unity when /a 
is put equal to unity. 

We shall next prove the following important proposition. 

If\ihe less than unity, and if (1 - 2/Mh. + h*)~* be expanded 
in a series proceeding hy ascending powers of h, the coefficient 
o/h* K?tK 6^Pj. 

Or, (1^2/^h + hT* = Pe+P,h + ... + P,h' + ... 

We shall prove this by shewing that, if H be written for 
(1 — 2fjJi + A*)~*, if will satisfy the differential equation 

For, since H= (1 - 2^ + A')"*, 



jp=l-2fih-i-h' 



Idff „. 

= -2/tff' + 3(l-/**)5'^ 
= -2/iH' + Sil-fi')hH\ 

= H*{l-2fth + V + h(ft-h)] 

=if'(i-M); 
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8 . ZONAL HAKMONICa 

= 3 (1 - ^A) fi* G* - h) -iiE\ 

3Mfi'+3{/x(l + /.*)-2/t»A}^ 

3/t{I?'-(l-2/iA + A.')fl»} 

= 0, since 1 - 2/iA + fe' == S^. 

Therefore | j(l - ;.^) f } + A |, (Afl) = 0. 

This may also be shewn as follows. 

If X, y, z be the co-ordinates of any point, / the distance 
of a fixed point, situated on the axis of z, from the origin, 
and R be the distance between these points, we know that. 



and that 



'•©=«• 



Now, transform these expressions to polar co-ordinates, 
by writing 

a? = rsin^cos^, y = rsintfsin^, j? = rcos5, 

and we get 

iJ' = r*-2/rcos^-f.«'«, 

and the differential equation becomes 



Digitized by VjOOQ IC 



ZONAL HARMONICS, 
or, putting cos ^ = /i, 
r 



/A, we see that 
= A' - 2/iA + 1 = 



Now, putting r = zh, we see that 
1 H 



R- 



/. the above equation becomes 

4. Having established this proposition, we may proceed 
as follows: 

Hpi be the coefficient of h* in the expansion of H^ 
jy= 1 +jp,A+jp^^ + ... +jp,fe* + ... 

Also, the coeflScient of h* in the expansion of 

i{<'-'*')3fl"|;{('-''')f- 

Hence equating to zero the coeflScient of h\ 
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10 ZaSLLL HABMONICS. 

Also p^ is a rational integral Tonctiou of /t. 
And, when /a = 1, ir= (1 - 2A + h^^ 

= l+A + /t' + ,..+A' + ... 

Or when /i = 1, p, = 1. 
Therefore p^ is what we have already denoted by P,. 

We have thus shewn that, if A be less than 1, 

If h be greater than 1, this series becomes divergent. 
But we may write 

since r is less than 1, 

Hence P, is also the coeflScient of A"<*+^> in the expan- 
sion of (1 — 2/xA + A*) "* in ascending powers of r when A is 

greater than 1. We may express this in a notation which is 
strictly continuous, by saying that 

Pi = P-(i+i)- 

This might have been anticipated, from the fact that the 
fundamental differential equation for P, is unaltered if 
— {% + 1) be written in. place of i ; for the only way in 
which i appears in that equation is in the coefficient of P,, 
which is i (i -f 1). Writing — (i + 1) in place of i, this be- 
comes — (i + 1) {— (i +!) + !} or (i + l)t, and is therefore 
unaltered. 
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SJONAL HARMONICS* H 

5. We shall next prove. that 

•where r* = a? + 3^ + a*. 

Let ^=(a^ + y* + «*)-*=/(^), 

and let k be any quantity less than n 

Then {3?+}f* + (z-ky}-^^/iz-k), 

and, developing by Taylor's Theorem, the coefficient of k' b 

(-')'i4^r«'t-»'r^^(;)- 

Also {x* + 2/» + («-;&)*}-* = (r'-2fe + A0"* 

since « = /^r, 

in the expansion of which, the coefBcient of k* is 

p. 

Equating these results, we get 

P^f^iy ^"^ 1(1] 

The value of P, might be calculated, either by expanding 
(1 — 2fih + h*y^ by the Binomial Theorem, or by effecting the 

differentiations in the expression (— 1)* ^ ^ •X*(~)» 

and in the result putting - = /it. Both these methods how- 
ever would be somewhat laborious ; we proceed therefore to 
investigate more convenieat expressions. 
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12 ZONAL HARMONICS. 

6. The first process shews, by the aid of Lagrange's 
Theorem, that 

Let y denote a quantity, such that 

h being less than 1, 

Theu 

1 
dy^ h, 1 

A y = /t + A(^^). 
Hence, by Lagrange's Theorem, 

therefore, differentiating with respect to /* and observing 
that 

+ 1.2...1VV 2 >/+'••• 
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ZONAL HA^MpNICa 13 

7. From this form of P, it may be readily shewn that 
the values of fi^ which satisfy the equation P^ = 0, are all real, 
and all lie between — 1 and 1. 

For the equation 

(ft*— 1)* = has t roots =1, and i roots = — 1, 

A -T- (jjf^Vf^Q hasi — 1 roots = 1, (t — 1) roots = — 1, and 
one root = 0, 

■5-5 (/A* — 1)' = has (t — 2) roots = 1, one root between 1 

and 0, one between and = — 1, and (i — 2) roots = —1, 
and so on. Hence it follows that 

T-, (jjf — 1)' = has s roots between 1 and 0, and - roots be- 
tween and — 1, if i be even, 

• 1 'I 

and — ^ roots between 1 and 0, —^ roots between and 

1, and one root = 0, if i be odd. 

It is hardly necessary to observe that the positive roots of 
each of these equations are severally equal in absolute mag- 
nitude to the negative roots. 

8. We may take this opportunity of introducing an im- 
portant theorem, due to Rodrigues, properly belonging to 
the Differential Calculus, but which is of great use in this 
subject. 

The theorem in question is as follows: 
Ifm he any integer less than i, 
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14 SX>^AL HAKMONICS. 

It may be proved in the following manner. 

If (aj*— 1)' be diflferentiated i -m times, tjien, since the 
equation 

has i roots each equal to 1, and t roots each equal = — 1, it 
follows that the equation 

has i — {i — m) roots (I. e. m) roots each = 1, and m roots 
each = — 1, in other words that (a? — 1 j** is a factor of 

We proceed to calculate the other factor. 

For this purpose consider the expression 

(x + a^) (x+a;) ... (x + a,) (a? + /3J {x + fi;} ... (^ + A)- 

Conceive this diflferentiated (I) i — m times, (II) i + m 
times. The two expressions thus obtained will consist of an 
equal number of terms, and to any term in (I) will corre- 
spond one term in (II), such that their product will be 
(cc + a^) {oo + a^) ... (x-\^a,) {x + l3^) (a? + /3,) ... (a? + y3,), i.e. the 
term m (II) is the product of all the factors omitted from 
the corresponding term in (I) and of those factors only. 
Two such terms way be said to be complementary to each 
other. 

Now, conceive a term in (II) the product of p factors of 
the form x + a, say x-\-a\ a; -h a" ... a? + a^^, and of q factors 
of the form x-\-^, say x + fi^, x-^ fi,,...x + I3^^y We must 
havej? + j = i — m. 

The complementary term in (I) will involve 

p factors x + 0, x + ff' ...x+ff^\ 

q factors x + a^, x + a^^ ... x +a(g). 
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ZONAL HAEMONICS. 15 

Now, every terra in (I) is of i -f m diraensions; We have 
accounted for p + j (or i — m) factors in the particular term 
we are considering. There remain therefore 2m fs^ctors to 
be accounted for. None of the letters 

a, a'\..o!^\ p^y ^/, •••A«)» 

can appear there. Hence the remaining factor must involve 
m a's and m ^'s, — say, 

There will be another term in (II) containing 

{x-¥^) (^ + /3") ... (aj + /3^') (a: + a>(^ + fltj ... (a? + aj. 

The corresponding term in (I) will be, as shewn above, 

{x-^ol) (a: + a") ... {x-¥c.^) ix + ^,) (oj + ^J ... {x + fij 

(a; + ,a) (aj + ,ot) ... (a; + «ot) (ar + ,^) (aj + ^) ... (aj-f^). 

Hence, the sum of these two terms of (I) divided by the 
sum of the complementary two terms of (II) is 

{x + ,a) (x + ,a) ... (x + ^ot) (aj-f,3) (x + ,$) ... {x+J5). 

Now, let each of the a'a be equal to 1, and each of the ^s 
equal to — 1, then this becomes {a^ — 1)'". The same factor 
enters into every such pair of the terms of (I). Hence 

Or ^^,^ = (a^ - 1)"* ^^ , to a numerical 

factor prh. 

The factor may easily be calculated, by considering that 
the coeflScient of x*^ in — 'jn^ — ^is2i(2i— l)...(i-fm+l), 

and that the coeflScient of a?*"*" in ^ri+^ — '^^ 

2i(2i-l).,.(i + m+l)(f + m)...(i-m + l). 
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16 ZONAX HARMONICS. 

Hence the factor is 

1 1.2... (i-m) 

(t + m) (i + m — 1) ... (i — m + 1)' 1.2... (t + m) * 

9. This theorem affords a direct proof that C-t-j (/i* — 1)', 

C being any constant, is a value of / (/a) which satisfies the 
equation 

from above, 
or 

|[(i-''')|.!|.('''-vl]+'('^^>{|(''-i)'}-«- 

Hence, the given differential equation is satisfied by put- 
ting/(M) = C7|^.(/**-l)'. 

Introducing the condition that P^ is that value of / (/tt) 
which is equal to 1, when /a = 1, we get 

-^'='2M.2...i5?^'~^)'- 

10. We shall now establish two very important proper^ 
ties of the function P< ; and apply them to obtain the develop- 
ment of P< in a series. 



Digitized by VjOOQ IC 



;50NAL HARMONICS. 17 

The properties in question axe as follows t 
Ifi and m he rmequal positive integers, 

And />^/* = 2?Ti' 

The following is a proof of the first property^r 

We have 

Multiplying the first of these equations by P^ the second 
by P^ subtracting and integrating, we get 

K{<^-''''f}*-M{<^-'**'f}'*' 

+ {i(i+l)-m(m + l)}jP,P„c?/i = 0. 

Hence, transforming the first two integrals by integration 
by parts, and remarking that 

i {i+ 1) — m (w + 1) = (i -m) (i+ m + 1)^ 
we get 

+ (i-m) (i+m + l)jP,PJlfi = 0, 



or 



since tlxe second term vanishes identically. 
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18 «ONAL HABMONICS. 

Hence, taking the integral between the limits — 1 and 
4- 1, we remark thc^t the factor 1 — /i' vanishes at both limits, 
and therefore, exc^t when i — m, or « + wi + 1 = 0, 

-1 
We may remark also that wo have in general 



/: 



i: 



*^ dpi. * dfJL 



a result which will be useful hereafter. 

11. We will now consider the cases in which 

i — m, or i 4- m + 1 = 0. 

We see that i+ m + 1 cannot be equal to 0, if * and m are 
both positive integers. Hence we need only discuss the 
case in which m = i. We may remark, however, that since 

P^ = P.(i+i), the determination of the value of i PJ d/ii will also 
give the value of I PiP^+i) dfL. 

The value of I P^dfi may be calculated as follows : 

(l-2M+A«j-i = Po-l-P,A + ...-fPA' + ...; 
.-. (l-Ji/iA + AT^ = (P,+ P,A+..,+P,/iV+...)V 
= Po* + P,«A»+...+P/A«+... 
+ 2P,P,A + 2P,P^«+ ... + 2P,P,A» + ... 

Integrate both sides with respect to /*; then since 
j{1^2pih + hYdfJL=^-^logil^2fih+hr), 
we get, taking this integral between the limits — 1 and + I, 

all the other terms' vanishing, by the theorem just pi o zed. 
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■ = r p;dii +h*f F*dft + ... + A" r p,w/t + ... 

Hence, equating coefficients of h", 



I 



'/■*-2T+T 



12. From the equation I P,P„, c?/a = 0, combined with 

• the fact^thatj.idien /i = l, P|==l, and that P^ is a rational 
integral function of /a, of the degree i, P^ may be expressed 
in a series by the following method. 

We may observe in the first t)lace' that, if m be any 

integer less than i, I fju^P^dfi = 0. 

■ • - • ^ 

For as P«, P^i... may all be expressed as rationarin- 

tegral functions of /a; of the degrees m, w — 1 ... respectively, 

it follows that /*"* will be a linear function of P^ and zonal 

harmonics of lower orders, /***"* of P^j and zonal harmonics of 

lower orders, and so on. Hence Ifi'^P/ifji will be the sum of 

a series of multiples of quantities of the form 1 PJPidii, 

m being less than i, and therefore I ^jJ^P^dfi* = 0, if m be any 
integer less than i. 
Again, since 

(l-2M + A')-* = -P. + -Px^+-+-P,^' + - 
it follows, writing — h for A, that 

(1 + 2M + A»r4 = P.-P,A+ ... + (- l)*P,;t*+ ... 

2—2 
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20 ZONAL HABM0NIC9. 

And writing — /a for /a in the first equation, 

Pq', P/...P/.., denoting the values which F^, P^.P^, 
respectively assume, when — /i is written for /*• Mence 
P^=zP^ or — P,, according as i is even or odd. That is, 
P, involves only odd, or only even, powers of <, according 
as i is odd or even*. 

Assume then 

P.^A.tif^'A^pf^-^... 

Our object is to determine -4.^, -4^.... 

Then, multiplying successively by /a*^, /a*^, ... and inte- 
grating from — 1 to +1^ we obtain the following system' of 
equations : 

A A A 

2i-1^.2i-3^-^2i-.25-l^-~"* 
2i-.3^2i-6^-^2i-.2«-3^- "' 



-^« J -4*=t_i -L •^<-y 



+ €>-• o« Q + »»' + o7 i:i =- + ..• — 0. 



2i-25-1^2i-2s-3^*'*^2i-45-l 



And lastly, since P^ = 1, when /* = 1, 

-4, + ^^ + ..• + -4^ + ..• - 1 ; 
the last terms of the first members of these several equa- 
tions being 

13. The mode of solving the class of systems of equa- 
tions to which this system belongs will be best seen by 
considering a particular example. 

* This is fdso eyident, from the fact that P| is a constant multiple of 
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Suppose then that we have 

a+a o+a c+a ' 

X y , « 1 

* -r 



a + a> 6 + 01 C + & 

From this system of equations we deduce the following, 
being any quantity whatever, 

_x_ ^ g _ 1 (g-g)(g-/3)(a+fl))(64-fl))(G+o)) 

For this expression is of — 1 dimension in a, b, c, a, /8, y, 
0, ft); it vanishes when = a, or = fi, and for no other 

finite value of 0, and it becomes = — , when d = o). 

We hence obtain 
_./ ^a^f y a. ^ \ 1 (g-a)(g-/3)((»+co)(6+ft))(c+ft>) 

and therefore, putting ^ = — a, 

1 (a + g)(a-f/3)(a + Q))(& + ft>)(c + co) 
^■"a> (a-.6)(a-c) (a>-a)(ft>-^) ' 
with similar values for y and z. 

And, if o> be infinitely great, in which case the last 
equation assumes the form a: + y + ^r = 1, we have 

_ (a + tt)(a + g) 

with similar values for y and z. 

14. Now consider the general system 

-?'— + ^i^^+. ..+-%-+. ..=0, 
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^' ■ '^'^ +...+ _^t^+...»o, 



o, + o^ o« + <»Mi. a^ + a.. 



+ ••• — f 



the number of equations, and therefore of letters of the 
forms a? and a, being — ^ if i be odd, ^ -f 1 if i be even ; and 

i — 1 . 

the number of letters of the form a being — ^ if i be odd, 

and ^ — 1 if I be even. 



We obtain, as before^ 



:+...+ 



Q> (G)-a,)(a>-a^J...((»-a^J...(a,+^)(a^,+^)...(a^jj.+ ^)... 
and, multiplying by a,_2,+ ^, and then putting ^ = — a^^, 

(a^ - a,) (a,.„ - a^ ... {a^^ - a^ or aj ' 

15. To apply this to the case of zonal harmonics, we see, 
by comparing the equations for x with the equations for A, 
that we must suppose « = qo ; and 

(^i = h a<lj = i — 2,...a,_2, = i'-25... 
a^ = i — 1, a,_a = » — 3,...a^j, = i — 2«— 1..- 
Hence 
. (2^-2^-l)(2^-2g-3)...{2(^-2j?)-l|.., 

*-«• "" (-2s){-(25-2)}...{(i-2^-l) or (i-2^)} 



= (-!)' 



.(2i-25-l)(2{-25-3)...{2(i-2»)-l|... 
2s(2«-2)...3x2,4..,Ci-2s.-l)or(i-25) 
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Or, generally, if i be odd, 

. ^ (2^•-l)(2^-3)...(^ + 2) 
* 2.4...(i-l) ' 

. ^_(2i--3)J2i-5)_* 

♦-• 2.4..,(i-3)x2 ' 

. _ (2;-5)(2f-7)...(t-2) 

*-• 2.4...(f-5)x2.4 ' 



^^"^ ^^ 2.4...(t-l)- 
And, if i be even, 

. _ (2t-l)(2t-3)...Ct+l) 

. (2t-3)(2V-5)...(t-l) 

*-• 2.4...(i-2)x2 ' 

. (2^^ - 5) (2i - 7) ...(t - 3) 

•-^ 2.4...(t-4)x2.4 ' 



_ i (»-l)(t-3)...l 

^•-^"•^^ 2.4...i 

We give the values of the several zonal harmonics, from 
P, to P^dnclusive, calculated by this formula, 

P. = l, 

P,= /*, 

_^ 3/t*-l 

~ 2 ' 
p 6 , 3 
■P. = 2/*-2^ 
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* 2.4'* 2. 2'*'^ 2.4 

* 8 ' 

p 9.7 , 7.5 ,^5.3 

_ 63/<.'-70/it*+15/* 

-: 8 ' 

P _ n.9.7 , 9.7.5 « 7.5.3 , 5.3.1 
' 2.4.6'* 2.4x:s'*"^2x2.4'* 2.4.0 

_ 231m* - 315/t* + lOS/i* - 5 
16"' ' 

p 13.11.9 , 11.9.7 , 9.7.5 , 7.5.3 
' 2.4.6'* 2.4x2'* ■''2x2.4'* 2.4.6'* 

_ 429/i" - 693/*' + 315/*' - 35/t 
16 

p _ 15.13.11.9 , 13.11.9.7 . 11.9.7.5 « 
• 2.4.6.8 '* 2.4.6x2'* ■^2.4x2.4'* 

9.7.5.3 , 7.5.3.1 
2x2.4.6'* "'■2.4.6.8 

_ 6435/a'' - 12012/t« -t- 6930/** - 1260/1,* + 35 
128 * 

p 17.15.13.11 . 15.13.11.9 .13.11.9.7 , 
•~ 2.4.6.8 '* 2.4.«x2 '*'^2.4x2.4'* 

11.9.7.5 . 9.7.5.3 , 
~2x2. 4. 6^* ■'■2.4.6.8''* 

_ 12155/1,* - 25740/i,^ + ISOlS/t' - 4620/t' + 315/t 
128 ' 

p _ 19.17.15.13.11 „ 17.15.13.11.9 , 15.13.11.9.7 , 
""2.4.6.8.10 '*■ 2.4.6.8x2 '*''*2. 4. 6x2. 4'* 
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13.11:9.7.6 ,. 11.9.7.6.3 , 9.7.5.3.1 



^ "^ O N^ Q A, CK q/* 'i 



2. 4x2. 4. e'^ '2x2.4.6.8'^ 2.4.6.8.10 

_ 46189/i^"- 109395/t^-i- 90090/i"- 30030/^"+ 3465/i»- 63 

256 

It will be observed that, when these fractions are reduced 
to their lowest terms, the denominators are in all cases 
powers of 2, the other factors being cancelled by correspond- 
ing factors in the numerator. The power of 2, in the 
denominator of P^ is that which enters as a factor into the 
continued product 1 . 2...i. 

16. We have seen that / /a"P,.(7/a = 0, if m be any 
integer less than i. 

It will easily be seen that if w + i be an odd number, the 
values of jfi'^P^.d/j, are the same, whether /a be put=sl or 

— 1; but if m+t be an even number, the values of l/A^P^.d/A 

corresponding to these limits are equal and opposite. Hence, 
{m + i being even) 

J%«P,.rf/i=2jVi'i.^M> 

and then I /A^P^.d/^ssO^ if m = t — 2, t — 4 

Jo 

We may proceed to investigate the value of / /**" P. . i/i, 

Jo 
if m have any other value^ For this purpose, resuming the 
notation of the equations of Art. 13, we see that, putting 
Ossm + l, and a> = 00 , we have 



a^ + m+l a,^ + m + l a^^ + m+l 

_ (m + 1 -- gQCm -f 1 - O ... (m + 1 - Oj^^) . 
{a^ + m + l){a^^ + m + 1) ...(a,^ + m+ 1) , 
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26^ ZONAL HABMOKICS; 

and tBerefore, patting Xf^^A^^.f a, » t . . ., a, = t' — 1 . . . , 
we get 
1 A A A 

o'*"-^'-'^'*"»+m+l"*"i-2+m+l ■•"•••+ t-28+m+l'^'- 
(m-i+2){m-i + 4.)...(m-l) if.t,.-jd 
(m + 1 + l)(m + i - 1) ... (»» + 4) (m + 2) " * "^ *^"' 
and (m-t + 2)(m-t+4)...m jfi-bccTcn 

In the particular case in which m = «^ we get 
£^'^'^^° C2U 1) (L-- 1) .^.^-1^4) (^ + 2) ^^'^^' 

*^^ ° (2^+ 1) (2»-lf.'.'.'^-+3) (» + !) (^'"""°)-' 

1 7. We may apply these formulaB to develope any positive 
integral power of /^ in a series of zonal harmonics, as we 
proceed to shew. 

Suppose that m is a positive integer, and that fi"* is de-- 
veloped in such a series, the coefficient of P, being C^, so 
that 

then, multiplying both sides of this equation by P, and inte- 
grating between the limits —1 and 1, all the terms on the 

right-hand side will disappear except I C^i^^dfju, which will 

2 ' 

become equal to q. .. C^. 



Hence <^* = ^^J A*"^«^/^' 



which is equal to 0, if m + 1 be odd. Hence no terms appear 
unless m + * be even. In this case we have 



J Q . . 
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Hence the formtila just investigated gives 
r (9!,!^ (w-f+2)(w»-t + 4)...(OT-l) 
^*~^^*'^^''(m+» + l)(m + «-l)...(m + 4)(m + 2> 

if i be odd, and 

^ - /9.- J. ,x (m-i+2) (m - i+ 4) .., m 

if « be even. 

Therefore if m be odd, 

«_,9 .,v 2.4.6...(ot-1) p . 

/* -t^»» + A;(2m+l)C2m-l)...(t» + 4)(m+2)^~+"* 

^ (m + 4)(m + 2)-^»^m + 2 >• 
If m be even, 

/t - ^^"^^ ^-'(2w+ 1) (2m-l) ... {m + 3) (»! + 1) -^^ "*" "* 

I 5 "^ P I ^ P 

Hence, putting for m successively 0, 1, 2 ... 10, we get 

/** = 5g-gP,+ gP,= gP, + gP„ 

'^ '7.5 '^S* 
5 »^6 *' 

^£p+*P+lp 



33*«^7*»^5 •' 
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"63 •"'■9 •■^7 " 
'* ^^''iS.ll.a.T '^^11.9.7 *^^ 9.7 »^7' 

-ii P+?fp^lOp4.1p 

"231 • 77 *'''21 » 7 "^ 

■ ^ 2.4.6 ' ,, 4.6 . 7_6_p . 3p 

'^ "*^ 15. 13. 11. 9^'^" 13. 11.9^'^ '11. g'' • + 9'^" 

"429 '^39 » 33 »^3 »' 

, H, 2.4.6.8 p 4.6.8 p 

'* ;^^ 17. 15. 13.11. 9'^«^^'* 15. 13.11.9^* 

■•■ * iOiTg ^^ "*■ ^ TO ^» ■*■ 9 '^'" 

128^ 64 48 40 1 

6435 « ^ 495 " ^ 143 "^^ ^ 99 » ^ 9 •' 

2.4.6.8 p ,,g 4.6.8 „ 
'* " ^ 19.17.15.13.11 •^17.15.13.11'^' 

I n ^'^ P t 7 ^ p I q ^ p 
■^ 15.13.11 »"^ '13,11 »"** 11 >' 

_ 128 192 16 56 3 
" 12155 • "^ 2431 ' "*■ 65 » "^ 143 » "^ 11 »' 

2.4.6.8.1 p y 4.6.8-.10 
'^ "21.19.17.16.13.11 "^19.17.15.13.11 » 

_6^8J0_ 8.10 p . 5_10_ 1 p 
+ ^** 17.15.13.11 ■^•+^15.13.11 '^'^^iS.ll'^'^ll"^' 

_ 256 p .128 p . 32 48 50 Ip 



46189 '•^2717 '^ 187 ,'^ 143 *^ 143 »^ 11 
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18. Any zonal harmonic P, may be expressed in a finite 
series of cosines of multiples of 0, these multiples being 
id, {i- 2)0.... Thus 

therefore, writing cos for /a, and observing that 

l-2cos^A+A« = (l-Ae^iO(l-*«"'^^0. 
we obtain 

(l_A€^^^-i(l-A€->^r^ = P, + P^A + ... + P;i' + .^ 

or 

whence, equating coefficients of h*, 
J, 1.3...(2i-l)_ ..^1.3...(2i-3)l- ,. -,. 

^1.3...(2t-5)1.3„ ,. ,,.^ 

the last tenn being ] o V ~ f if* ^ even, and 
1.3...(i + l) 1.3...(i-2)^ „..., ,. 

19. . Let us next proceed to investigate the value of 

/ P| cos m0 sin ^ d0. 
•'0 
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This might loe done, by direct integration, from the above 
expression. Or we may proceed as follows. 

The above value of P^ when multiplied by cos md sin 
(that is by 2 t^^^ (^+1) ^ — sin (m — 1) 0]) will consist of a 

series of sines of angles of the form {% — 2n ± (m ± 1)} 0, that 
is of even or odd multiples of ^, as t + m is odd or even. 
Therefore, when integrated between the limits and tt it 
will vanish, if i + mhe odd. We may therefore limit our- 
selves to the case in which i + m is even. 

Again, since cos m0 can be expressed in a series of powers 
of cos 0, and the highest power involved in such an expression 
is cos""^, it follows that the highest zonal harmonic in the 

development of cos m0 will be P^. Hence i P^ cos m0 sin d0 

will be = 0, if m be less than t. 

Now, writing 

P,= Cicos;^+(7^cos(i-2)^ + ..* 

we see that P^ cos m0 sin d0 will consist of a series of sines 
of aiigles of the forms (m + 1 + 1) 0^ (m + 1 — 1) ^ ... down to 
(m — t — 1)^, there being no term involving m^, since the 
coefficient of such a term must be zero. Hence 



/• 

Jo 



P,cos«i^sin^(?^, 

'0 

will consist of a series of fractions whose denominators in- 
volve the factors m + i + l, w + i — 1...77^ — t — 1 respectively. 
Therefore when reduced to a common denominator, the result 
will involve in its xienominator the factor 

(m + t + l)(m + t-l)...(m + l)(m-l)... (m-t-1) 
if m be even, and 

(m + * + l)(w + *-l)...(m + 2)(m-2)...(m-t-l) 
if m be odd. 

For the numerator we may observe that since 
P^ cos m sin 0dd 



■/: 
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•vanishes i£ m he less than t, it must ii^volve thU 6»ctors 
w — (i — 2), w— (i — 4)... »n + (t — 2), and that it does nOt 
change sign with m* Hence it will involve the factor 

{m-(t-2)}{m-(i-4)}.,. (m-2)i»«(w + 2)...(m + t-2) 

if w T^e even, and ... 

{m - (t - 2)} {m- (t - 4)} ... (w - 1) (m + 1) ,., (m + » - 2) 

if m be odd. 

To determine the factor independent of m, we may pro- 
ceed as follows : 

P, = Cicos 10 + C^ cos (t- 2) ^ + ... ; 
/. P,cosm5 = 5 C^{cos(m + i)^ + cos (m — t) 6] 

+ 2 CL, {cos (w + 1- 2) ^ + cos (m- i + 2) ^} + ... ; 

.-. P, cos m^ sin O^-rC^ {sin (m + i + 1) ^-sin (w + i- 1) ^ 



+ sin (m — «+l)^ — 8in(m-f— 1)5} 
+ ict,{sin(m + t-l)5-sin(w + t-3)5 

+ sin(w-t + 3) 5- sin (w - iH- 1) 5} + ... ; 
/. I PjCOS w5sin5(?5 

Ifif 1 - 1 + 1. - M 

2[wi + i+l m + t— 1 m — 1 + 1 m — t — lj 

• +^.|^ L_+ \ i— 1 + ... 

2 \m + i — 1 ?^^ + i — 3 wi-t-f3 m — i+lj 
■•" ^'^. 1" m'- (i - 1)' ^^m'- {i-m"" '" 
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Now, when m is very large as compared with t, this be^ 
comes 

since C, + (7^, + ••• = 1, as may be seen by putting = 0. 

Hence / P. cos m^ sin ^ d^ tends to the limit -^asm 

Jo ^ . ^ 

is indefinitely increased. 

The value of the factor involving m has been shewn 
above to be 

{m-(i-2)}{m-(t-4)}..,(m-2)m»(m + 2)... (m + i-2) 
{m - (t + 1)} {m - (i - 1)} ... (m - 1) (m + 1) ... (m + 1 + 1) 

if m be even, and 
{m-(t-2)}{m-(i-4)}...(m-l)(m + l)... (m + t-2) 
jw - (i+ 1)} {m - (t - 1)] ... (m - 2) {m + 2) ... (m + 1 + 1) 

if m be odd. 

Each of these factors contains in its numerator two factors 
less than in its denominator. It approaches, therefore, when 

m is indefinitely increased, to the value -j • Hence 



/. 



r 

P^ cos md sin dO 



{m-(t-2)}{m-(i-4)]...(m-2)m»(m+2)...{m+(t-2)} 
- '^{w-(i+l)}{m-(i-l)}...(m-l)(m + l)...tm+(t + l)} 
if m and i be even, and 

[m-(i-2)}{m-(t^4)}...(m-l)(m+l)...{m+(t-2)} 
^{w-(;+l)}{i^-(4-l)}..,(m-2)(m+2)...{m+(t+l)| 

if m and « be odd* 

In each of these expressions % may be any integer such 
that m — « is even, i bemg not greater than m. Hence they 
will always be negative, eiscept when % is equal, to m. 
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20. We may apply these expressions to develop cos mO 
in a series of zonal harmonics. • 

Assume 

Multiply by P« sin 0^ and integrate between the limits 
and IT, and we get 

g {m-(i-2)}{»t-(t-4)}...{m+(t-2)} 2 
{m-{i + l)}[m-(S-l)} ...t»» + (« + l)} 2» + l *• 

Hence 

B- (-i I i J^-(^-2)}{m~(/-4)}...{m + (i-2)} 

Hence, putting m successively = 0, 1, 2, ... 10, 
cosO^ = P^; 
cos^ = Pj; 

3 * 3 •' 



3.5 



-l;j^^. 



"35 « 21 * 15 •' 
F. H. 
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r/i -i-i 2.4.6.8 „ ^ 4.6 „ 



-^sV' 



-128p_8 1 

~ 63 » 9 » 7 " 

«.«fifl- i« 2. 4. 6'. 8. 10 p 

**'^^ — ^^-1.1. 3. 5. 7. 9.11. IS-^* 



4.6V8 6' p 1 

''1.3.6.7.9.11 « 3.5.7.9 « 5.7 * 
_512 384 £p_J. p. 
231 • 385 * 21 * 35 •' 

^a ir 2.4.6.8.10.12 j, 

co87^=-16_j^ 3 g g ^^ J3 jgP, 

4.6.8.10 p ^ 6.8 p 3^ p- 

1.3.5.9.11.13 » '3.5.9.11 » 5.9^ 
_1024 128 112 J^ 
"429^ 117 » 495 » 15 »' 

o^ ,» 2.4.6.8M0.12.14 „ 

cos8d=-17_j J 3 5 7 9 jj j3 j5 ^yP, 

4.6.8M0.12 6.8M0 p 

1.3.5.7.9.11.13.15 • 3.5.7.9.11.13^ 

K 8* p__L_p 

^67779711 » 7.9 • 

_ 16384 p _4096p_J56 p 64 1 p 

~ 6435 * 3465 • 1001 * 693 « 63 "' 
Qfl- ,r> 2.4.6.8.10.12.14.16 
cosyp ^^-1.1.3.5.7.11.13.15.17.19 • 

4.6.8.10.12.14 p 6.8.10.12 

^^1.3.5.7.11.13.15.17 ' 3.5.7.11.13.16^» 

_ _8^10_ _1_ 

'5.7.11.13 » 7.11 ' 
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_ 82768 3072 p 128 „ _i6^p _ J^ p . 
"12166 • 2431^^455 • 143 » 77 »' 
in/) _ 91 2. 4.6.8.10'. 12. 14. 16. 18 p 

coslO(/ ^^_ 1.1. 3. 6. 7. 9. 11. 13. 16. 17. 19. 21 »• 
y 4. 6. 8.10'. 12. 14. 16 p 
^^1.3.5.7.9.11.13.16.17.19 ' 
- 6.8. 10*. 12.14 p 8.10*. 12 p 

"3. 5.7. 9. 11. 13. 16. 17 • ^^6. 7. 9. 11. 13. 15 * 

5 10* p l_p 

7.9.11.13 » 9.11 • 

131072 p 32768 p 512 p 128 p 500 p 
~ 46189 » 24453 ' 1683 • " 1001 * 9009 ' 

99 • 

21. The present will be a convenient opportunity for 
investigating the development of sin d in a series of zonal 
harmonics. Since sin ^ = (1 — /t*)i it will be seen that the 
series must be infinite, and that no zonal harmonic of an odd 
order can enter. Assume then 

sind=i 0,P,+ C,P, + ... + CiP,+ ... 

i being any even integer. 

Multiplying by P«, and integrating with respect to ft 
between the limits — 1 and + 1, we get 

f'p,sin^(?/t = 2^^Ci; 

... Ct=^^i±lp P^anOdfi 

supposing P( expressed in terms of the cosines of and its 
multiples 

3—2 
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Hence, putting i = 0, 

1 3 

Putting i = 2, and observitig that P, = t + t cos 25, 



^•-ii: 



(1 + 3 cos 20) (1 - cos 20) ,^ 
(j^ff 

4 



= ^r|l+2co8 25-|(l + cos4^|rfd 

For values of ^ exceeding 2, we observe, that if we write 
for Pi the expression investigated in Art. 18, the only part 

of the expression I P* (1 — cos 20) d0 which does not vanish 

will arise either from the terms in P< which involve cos 20, or 
from those which are independent of 0, We have therefore 

_ 2t + l 1.3. ..(«-!) 1.3. ..(»•- 3) 
4 '2.4... t 2.4... (»-2) 

[' /iri + *±1 2 cos 2tf) (1 - cos 2^ d^ 

2t + 11.3...(t-l)1.3... (t-3) fi-1 i+l\ 
"^ 4 2.4... » 2.4...(t-2)'^l « » + 2/ 
2t + 11.3... (t-1) 1.3... (»-3) 
'^ 2 2.4...«(i + 2)2.4...(»-2)»" 

Hence sm^ = Jp,-g|p,-... 

(2t + l)Trl.3... (i-1) 1.3... (/-3) ^ 

2 2.4...«(» + 2)2.4... (t-2)*-^*""' 
i being any even integer. 
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dP 

22. It will be seen that -7-* , being a rational and integral 

function of fi*'\ M*"^..., must be expressible in terms of 
-fLi> -fU««» To determine this expression, assume 

dP, 
dfi 



■ = CL,P^, + CL.P^3 + - + CLP^+,.. 



then multiplying by P^, and integrating with respect to fi 
from — 1 to + 1, 

Now, since I >m, 

■■■fj:^di,-[p,Pi'-iPj>r-% 

since either m or i must be odd, and therefore either P„ or 
P, = — 1, when /t = — 1 ; 



_2_ 
2m+l 



^ = szrr-^G^or CL=2m + l; 



.•.^'=(2i-l)P^, + (2i-5)P... + (2i-9)P^ + ... 

23. From this equation we deduce 

Pi~P^«-(2»-l)//i.i^/*> 

the limits /* and 1 being taken, in order that P, - P,.^ may 
be equal to at the superior limit. 
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Now, recurring to the fundamental equation for a zonal 
harmonic, we see that 






dP^, 



24. We have already seen that | PtP^ dj^^O, % and m 

•'-1 .... 

being different positive integers. Suppose now that it is 
required to find the value of J P,P«»^/a. 

We have already seen (Art. 10) that 

j^P.P«rf/*= (i_«,)(^ + m+l) • 

And, from above, 

r'^-r* I 1 fOT(m+l)p/p _p ^ 
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25. We will next proceed to give two modes of ex- 
pressing Zonal Harmonics, by means of Definite Integrals. 
The two fexpressions are as foUows : 

* ttJo {/*±(^"-l)*cos^l**^' 

P,=l|'{/,±(^«-l)Jcos^yrf^. 

These we proceed to establish. 
Consider the equation 

d^ 1 



-r- 



a-ftcosSr (a^^i^y 

The only limitation upon the quantities denoted by a 
and h in this equation is that 6* should not be greater than 
a\ For, if 6" be not greater than a* cos S^ cannot become 

equal to t while ^ increases from to tt, and therefore the 

expression under the integral sign cannot become infinite. 

Supposing then that we write z for a, and V— 1 p for 6, 
we get 

1 r* d^ _ 1 

'irJo ^-V^/DCOSS:" («'+/}')*' 

We may remark, in passing, that 

Jo « — V—l/OCOSSr ^0 Z + ^ — lpCOH^ 

zcF^ 



Jo i 



Jo «• + /)* COS^Sr' 

and is therefore wholly real 

Supposing that p*=^a? + jf, and that aj'+y* + a' = r*, we 
thus obtain 



1 f' ^ 



d^ 1 



TTJo « — V-lpCOSSr ^ 
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Differentiate i times with respect to z, and there results 
1 cf [• ^ (f 1 P, 



d'f' ^ ^^l^-Py lyi 2 3 i 

/OCOSSr 



IT 

TT p ^^ (-1)' <p r ^ 

Hence ■?. = — j" 273^:?^? Jo^-V^ 



TT Jo (a-V~l 



(a-V-lpcosS-r'* 
In this, "write iw for «, and (1 — /t**)* r for p, and we get 

which, writing 9r — ^ for Sr, gives 
IT* d!^ 

26. Agaiu, we have 

1 1 f' ^ d^^ 

(a*— J*)* TT Jo Ct — 6cOSi|r' 

In this write 1 -/^Afor a, and ± (/a'-1)* fefor6,which is 
admissible for all values of h from up to /*— Ott* — 1)^ and 
we obtain, since a* — 6* becomes 1 — 2/A + A*, 

1 1 r» di^r 

(l-2M + A")*""^Jo l-M+OiA*-l)*Acosi|r 

^if' ^f: . 

^rjo 1-(/ea±(/ia'-1)*cos^}A' 

=-^(^^[l+{At±OA«-l)*cost}A+... 

+ {itt ± (^'-l)* cos ^}'A' +...]. 
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Hence, equating coefficients of h\ 

Pi^W {/* ± (ifi*-l)*C08^}*rf^. 

The equality of the two expressions thus obtained for P, is 
in harmony with the fact to which attention has already 
been directed, that the value of P, is unaltered if — (i + 1) 
be written for i. 

27. The equality of the two definite integrals which 
thus present themselves may be illustrated by the following 
geometrical considerations. 




Let be the centre of a circle, radius a, C any point 
within the circle, PCQ any chord drawn through (7, and let 
00=6, GOP^%GOQ = ^. Then aP» = a* + 6*-2a6cos^, 
G(^ = a* + 6^ — 2a& cos '^. Hence 

(a'^ + 6'-2a6cos^)(a* + 6*-2aScosi^)=:(a*-6y; 

sin^c?^ sinifrrfifr _ 

•':a« + 6>'-.2a6cos^"*"a*+6«-2a6cosi|r"" ' 
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42 ZONAL HARMONICS. 

Again, since the angles OPC, OQC are equal to one 
another, 

BixiL0_ sm OPC _ sin OQC sin^ 
CF OU ^ 00 " CQ ' 

. sin^ sinifr 

"(a'+6"-2a6cos^)i"'(a* + 6'-2a6co8^)4' 

whence - ^ it—^o 

(a"+6*-2a6cos^)4"^(a' + 6*-2a6cos^)4 

In this, write a^ + V^fi, 2a& = ? (ft*— 1)*, which gives 
a^ — 5« r= 1^ and we get 

We also see, by reference to the figure, that as ^ in- 
creases from to 9r, -^ diminishes from 9r to 0. Hence 

/: i.t(.--t'co.^r -/>*»'-'^'""^''''^- 

28. From the last definite integral, we may obtain an ex- 
pansion of P| in terms of cos and sin 0. Putting /a = cos 0, 
we get 

P| = 2^J^[{cos + V^ cos yft sin 0)]* 

+ {cos - V^ cos -^ sin ^j*] d-^ 

«- n(cos^'-^^^cos"t(cos^*-*Csin^)"+... 
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+ (_!)« i{i-l).'-(i-^m + l) ^^ ^^„ ^^^ ^^^ ^^ ^^ 

. , %{%-!)...{%- 2m + 1) (2to- 1) (2m-3)...l 
1.2.. .2m 2m(2m-2)-..2 

_ f(t-l)...(t-2m + l) 
(2.4...2m)» ' 

.-. P,= (cos ff)*-^-^^^ (cos ^" (sin ^»+ ... 

4.(-ir ^'(^'-g-^;-;;^^^) (cosgr''(sing)'"+... 
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CHAPTER III. 

APPLICATION OF ZONAL HARMONICS TO THE THEORY OF 
ATTRACTrON. REPRESENTATION OF DISCONTINUOUS 

FUNCTIONS BY SERIES OF ZONAL HARMONICS, 

1. We shall, in this chapter, give some applications of 
Zonal Hai-monics to the determination of the potential of a 
solid of revolution, symmetrical about an axis. When the 
value of this potential, at every point of the axis, is known, 
we can obtain, by means of these functions, an expression 
for the potential at any point which can be reached from 
the axis without passing through the attracting mass. 

The simplest case of this kind is that in which the 
attracting mass is an uniform circular wire, of indefinitely 
small transverse section. 

Let c be the radius of such a wire, p its density, k its 
transverse section. Then its mass, M, will be equal to 2'rrpck, 
and if its centre be taken as the origin, its potential at any 

point of its axis, distant z from its centre, will be j-. 

Now, this expression may be developed into either of the 
following series : 

zT 27^2Az* -^^ ^^ 2A...2i ^« + -}-/A 

We must employ the series (1) if be less than c, or 
if the attracted point lie within the sphere of which the ring 
is a great circle, and the series (2) if z be greater than c, 
or if the attracted point lie without this sphere. 
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Now, take any point whose distance from the centre is r, 
and let the inclination of this distance to the axis of the 
ring be 0. In accordance with the notation already em- 
ployed, let cos O^/JL. Then, the potential at this point will 
be given by one of the following series : 

^^ ^ 2.4.0... 21 -^^^c'* V**^ ^' 

^^ ^ 2.4.6...2i ^«'^ + ---|--^^^ 

For each of these expressions, when substituted for V, 
satisfies the equation V*F = 0, and they become respectively 
equal to (1) and (2) when is put = 0, and consequently 
r = z» The expression (2') also vanishes when r is infinitely 
great, and must therefore be employed for values of r greater 
than c, while (1') becomes equal to (2') when r = c, and will 
therefore denote the required potential for all values of r 
less than c. 

These expressions may be reduced to other forms by 
means of the expressions investigated in Chap. 2, Art. 25, viz. 

or P,= i fOi + J^^iF^ cos ^|ry*''d^{r. 

Substitute the first of these in (1') and (observing that 
fir = z) we see that it assumes the form 

M f' ( l {^+(g^-y^)*cos^}' 
J 1 2 c' 



TTC 
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which is equivalent to 



Mr 



(;5 + (-5««r')4cos^H4' 



The substitution of the last fonn of P^ in the series (2') 
brings it into the fonn 

Mr ( 1 1 c' 

1.3 ^ \^ 

■*■ 2 . 4 {« + (i5» - r*)i cos^}* " • • -J "^^ 

which is equivalent to 



Mr 

TT Jo 



[{« + (^-/^)^cos^}« + c^i* 

2. Suppose next that the attracting mass is a hollow shell 
of uniform density, whose exterior and interior bounding 
surfaces are both surfaces of revolution, their common axis 
being the axis of z. Let the origin be taken within the 
interior bounding surface ; and suppose the potential, at any 
point of the axis within this surface, to be 

Then the potential at any p«int lying within the inner 
boimding surface will be 

A,P, + A,P,r + AJP,r^ + ...+A,P/ + ... 

For this expression, when substituted for F, satisfies the 
equation V'F=0; it also agrees with the given value of 
the potential for every point of the axis, Ijring within the 
inner bounding surface, and does not become infinite at any 
point within that surface. 

Again, suppose the potential at any point of the axis 
without the outer bounding surface to be 
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Then the potential at any point lying without the outer 
bounding surface will be 

r r* T* r^* 

For this expression, when substituted for V, satisfies the 
equation V*F= ; it also agrees with the given value of the 
potential for every point of the axis, lying without the outer 
bounding surface, and it does not become infinite at any 
point within that surface. 

By the introduction of the expressions for zonal har- 
monics in the form of definite integrals, it will be found that 
if the value of either of these potentials for any point in the 
axis be denoted by ^ («), the corresponding value for any 
other point, which can be reached without passing through 
any portion of the attracting mass, will be 



^r^{« + («*-r^)*cos^}d^. 



3. We may next shew how to obtain, in terms of a series of 
zonal harmonics, an expression for the solid angle subtended 
by a circle at any point. We must first prove the following 
theorem. 

The solid angle, svhtended hy a closed plane curve at any 
pointy is proportional to the component attraction perpendicular 
to the plane of the curve, exercised upon the point hy a lamina, 
of unciform density and thickness, hounded hy the closed plane 
curve. 

For, if dS be any element of such a lamina, r its distance 
from the attracted point, 6 the inclination of r to the line 
perpendicular to the plane of the lamina, the elementary 
solid angle subtended by d8 at the point will be 

dS cos 
r" • 
And the component attraction of the element of the 
lamina corresponding to d8 in the direction perpendicular 
to its plane will be 

pk -^cos0, 
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p being the density of the lamina, h its thickness. Hence, 
for this element, the component attraction is to the solid 
angle as pA; to 1, and the same relation holding for every 
element of the lamina, we see that the component attraction 
of the whole lamina is to the solid angle subtended by the 
whole curve as pk to 1. 

Now, if the plane of ay be taken parallel to the plane 
of the lamina, and V be the potential of the lamina, its 
component attraction perpendicular to its plane will be 

— ^. Now since Fis a potential we have V'F=0, whence 

— V*F=0, or ^(•^)=0. Hence -r- is itself a potential, 

and satisfies all the analytical conditions to which a potential 
is subject. It foUows that, if the solid angle subtended by 
a closed plane curve at any point (a?, y, z) be denoted by 
a>, 0) will be a function of a;, y, z, satisfying the equation 
v'a) = 0. Hence, if the closed plane curve be a circle it 
follows that the magnitude of the solid angle which it sub- 
tends at any point may be obtained by first determining 
the solid angle which it subtends at any point of a line 
drawn through its centre perpendicular to its plane, and 
then deducing the general expression by the employment 
of zonal harmonics. 

Now let be the centre of the circle, Q any point on the 
line drawn through perpendicular to the plane of the 
circle, E any point in the circumference of the circle. With 
centre Q, and radius QO^ describe a circle, cutting QE in L, 
From L draw LN, perpendicular io QO. 

IjetOE^c,OQ==z. 

ThenJFi = (c*+^)4-;^, Oi«^«^^rf^{(c» + ^)»--^} 



I 



■e- 



(c* + «»)*■ 
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And the solid angle subtended by the. circle at 
A ON 



= 27r(l '—y 



To obtain the general expression for the solid angle sub- 
tended at any point, distant r from the centre, we first 
develope this expression in a converging series, proceeding 
by powers of z. This will be 

if z be less than c, and 

„ (lc« 1.3 c* , T,.1.3...(2i-l)c»^ ] 

if z be greater than c* 

Hence, by similar reasoning to that already employed, 
we get, for the solid angle subtended at a point distant r 
from the centre, 

^^X" e +2 c» 2.4 c* ^••* 

^ -^ 2.4...2i c"*"-' ■*"•••) 
if r be less than c, and 

flP.c' 1.3P/ . 1.3...(2t-l) P^,c" ) 

^'^ 12 "^ -IT* 7^+ - " ^~^' 2.4...2i r«"+-; 
if r be greater than c. 

4. We may deduce from this, expressions for the potential 
of a circular lamina, of uniform thickness and density, at au 
external point. For we see that, if F be the potential of 
such a lamina, h its thicfcn^ss, and p its density, we have for 
a point on the axis, 



-f=/,^.2.{l-^| 



F. H. 
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whence F«27rpfc{(c* +«')*-«} 

if M be the mass of the lamina. 

Expanding this in a converging series, we get 

if a be less than c, and 

.afflc' 1.1c* 1.1.3 c' 
''^~c"l2a 2. 4 a' 2.4.6 a' '" 

1.1.3...(2t-.3) c" 1 

~^"^^ 2.4.6...2i a""^--"] 

if a be greater than c. 

Hence we obtain the following expressions for the po- 
tential of an uniform circular lamina at a point distant r 
from the centre of the lamina : 

•^ c'l" -^'^^2 c 2.4 c» ^ 

,. 1.1.8...(2t- 3) P^r" •) 

" ^~ ^^ 2.4.6...2» c»-' "^ • ■ "I 

if r bo less than c, and 

^_ if fl Po^ _ 1^1 P,c« 1JL3 ^ _ 
'^~ c't2 r 2.4 ♦* "'■2.4.6 i* '" 

. ,,. 1.1.3...(2t--3) P^c" . ■» 
"^"^^ 2.4.6...2t, r«-^ +--| 
if r bo greater than c. 
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It inay be shewn that the solid angle may be expressed 
in the form 

27r-2 r ^ + (^'-O^^Q^^ ^0 
J [c*+ {^ + («« - r^)i cos ^l^ji ' 

and the potential of the lamina in the form 

5. As another example, let it be required to determine 
the potential of a solid sphere, whose density varies inversely 
as the fifth power of the distance from a given external point 
at any point of its mass. 

It is proved by tbe method of inversion (see Thomson 

and Tait's Natural Philosophy, Vol. 1, Art. 518) that the 

M 
potential at any external point P' will be equal to ,p, , 0' 

being the image of in the surface of the sphere, and M 
the mass of the sphere. We shall avail ourselves of this 
result to determine the potential at a given internal point. 

Let G be the centre of the sphere, the given external 
point. Join CO, and let it cut the surface of the sphere in A, 
and in GA take a point (7, such that CO . GO = GA^. Then 
O is the image of 0. 

Let P be any point in the body of the sphere, then we 
wish to find the potential of the sphere at P. 

Take as pole, and DC as prime radius, let OP — r, 
POG=e. AlsoletC^ = a, (70 = c. 

Let the density of the sphere at its centre be p, then its 

density aft P will be p -5 . Hence 

M^^irlLjy&mddrde, 

4—2 
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the limits of r being the two values of r which satisfy the 
equation of the surface of the sphere, viz. 

»^ + c'-2crcos^ = a», 

and those of being and sin"* - . 

c 

Hence, if r^, r, be the two limiting values of r, we have 



AT 1 .1 2ccos^ 

Also - + - = -^ s- 1 

r, r^ c —or 

1 1 



rjr. (? — a' 



'11 






»{c'cos»d-(c»-o'))i 



_, fg* -c* sin* g)^ 



.-. M = — '^ - -J i . -^ J / cos sin (a*- c* sm* ^) » 



dd 



a 

-1- 
c 



= -^"^ ^'^' ' cos ^ sin e (a» - c' sin» ff}^ dO 



_ 47r/3C /"si 

Now, if F be the potential at P, we have (see Chap. i. 
Art. 1) 

^d^{rV)_^ 1 d (. .dV\ 47rpc' 
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This is satisfied by F= — ji -^ . 

Assume then, as the complete solution of the equation, 

+ (a/+^)p, + .... 

It remains to determine the coefficients -4^, A^...A^,,.B^^, 

B^...Bi, so that this expression may not become infinite for 

any value of r corresponding to a point within the sphere, 

and that at any point P on the surface of the sphere it may 

M 
be equal to jy-p, where (XP : OP :: a : c, and therefore, at 

the surface, 

«.__ Mc 1 __4s irpc^c? 
ir~OP''S{G''-a:yr' 

And, at the surface, we have 
r»-2crcos^ + c'-a' = 0; 

1 _ 1 /1 2c cos g ^ 

. 4 Tpc'a' _2 ^pc- /P. PA / BA 

+ (Aj'+-ij Pj + ... identically. 

/ 4^rp^ \ 
V 3c"-o''^ V " 
and P,, P„...P,...i„ ^j..,^, aE = 0. 

Digitized by VjOOQIC 



54 AFPUCATIOX OF ZONAL EARUONICS 

Hence since P, = 1, 



•^•"Sc'-o'Kc'-a") I 



whence we obtain, as the expression for the potential at any 
internal point, 

6. We shall next proceed to establish the proposition that 
if the density of a spherical shell, of indefinitely small thick- 
ness, be a zonal sfiirface harmonic, its potential at any internal 
point will be proportional to the corresponding solid har- 
monic of positive degree, and its potential at any external 
point will be proportional to the corresponding solid harmonic 
of negative degree. 

Take the centre of the sphere as origin, and the axis of 
the system of zonal harmonics as the axis of z. Let b be the 
radius of the sphere, Si its thickness, U its volume, so that 
U= 4mVSb. Let CPi be the density of the sphere, P, being 
the zonal surface harmonic of the degree i, and C any con- 
stant. 

Draw two planes cutting the sphere perpendicular to the 
axis of z, at distances from the centre equal to f, f + df 
respectively. The volume of the strip of the sphere inter- 
cepted between these planes will be ^ U, and its mass will be 

Now ?= J/A, hence d^=bdiJL, and this mass becomes 
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Hence the potential of this strip at a point on the axis of Zy 
distant z from the centre, will be 

OTJ P, . 

which may be expanded into 

and ^5(p^ + p^^+... + P.^;+....),i^if;,>S; 

To obtain the potential of the whole shell, we must inte- 
grate these expressions with respect to /i between the limits 
— 1 and + 1. Hence by the fundamental property of Zonal 
Harmonics, proved in Chap. II. Art. 10, we get for the po- 
tential of the whole shell 

CU z* 
IT' — r W^ ^* ^^ internal point, 



CU_V 

2i + 1 z' 



^ at an external point. 



From these expressions for the potential at a point on 
the axis we deduce, by the method of Art. 1 of the present 
Chapter, the following expressions for the potential at any 
point whatever : 

CU Pr* 

Fi = ;t-. — = TTPi ^* a^ internal point, 

Vr CU Ph' , . 1 • . 

y* = -ET' — T -iir at an external pomt. 
2 2i + 1 r "■ ^ 

From hence we dieduce the following expressions for the 
normal component of the attraction on the point. 

Normal component of the attraction on an internal point, 
measured towards the centre of the sphere, 

dr WPL^*^ V^ • 
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Normal component of the attraction on an external point, 
measured towards the sphere, 

In the immediate neighbourhood of the sphere, where r is 
indefinitely nearly equal to J, these normal component at- 
tractions become respectively 



2i + l^ V" 2i + l b" 
and tbeir difference is therefore 

And writing for U its value, 4i7rh^B'b, this expression be- 
comes 

47rS6.aP,. 

Or, the density may be obtained by dividing the alge- 
braic sura of the noraial component attractions on two points, 
one external and the other internal, indefinitely near the 
sphere, and situated on the same normal, by iir x thickness 
of the shell. 

7. It follows from this that if the density of a spherical 
shell be expressed by the series 

G,P,+ C,P,+ CJ', + ... + C,P, + ..., 

Cq, Cj, Cg ... C, ... being any constants, its potential (FJ at 
an internal point will be'' 

^V b 3 b* +5 b' +- + 2i + l b'*' ^•"■J 
and its potential ( F,) at an external point will be 

In the last two Articles, by the word '^density" is meant 
"volume density," i.e. the mass of an indefinitely small 
element of the attracting sphere, divided by the volume of 
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the same element. The product of the volume density of 
any element of the shell, into the thickness of the shell in 
the neighbourhood of that element, is called "surface den- 
sity.** We see from the above that, if the surface density 
be expressed by the series 

the potentials at an internal and an external point will seve- 
rally be 

\ r 3 r^ 6 r 2i + 1 r^ / 

This variation in surface density may be obtained either 
by combining a variable volume density with an uniform 
thickness, as we have supposed, or by combining a variable 
thickness with a uniform volume density, or by varying both 
thickness and density. 

8. We have seen, in Chap. Ii., that any positive integral 
power of fi, and therefore of course any rational integral 
function of fi, may be expressed by a finite series of zonal 
harmonics. It follows, therefore, that we can determine the 
potential of a spherical shell, whose density is any rational 
integral function of /x-. 

Suppose, for instance, we have a shell whose density 
varies as the square of the distance from a diametral plane. 
Taking this plane as that of a?y, the density may be ex- 

pressed by pfi\ or p^. We have seen (Chap. n. Art. 20) 

that 

,*' = |(1 + 2P0. 

Hence, by the result of the last Article, the potential 
will be 

P "^ VS "^ ^ "^) ^* ^^ internal point, 
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p o^ ( ~ + ¥ ~|~ J at an external point ; 

ot, since P/' = -^ — r =» — er^» ^^ obtain 
p-^ (s + H ^ — ) ^^^ *^® potential at an internal point, 
Po 1 "^Tf—^"*" -rjf for that at an external point. 

9. As an example of the case in which the density is re- 
presented by an infinite series of zonal harmonics, suppose we 
wish to investigate the potential of a spherical shell, whose 
density varies as the distance from a diameter. Taking this 
diameter as the axis of z, the density will be represented 
by p sin 0, or p (1 —/x-*)*. We have investigated in Chap. II. 
Art. 21, the expansion of sin in an infinite series of zonal 
harmonics. Employing this expansion, we shall obtain for 
the potential 

ZeEilp-Lp'l 1.3...(t-l) 1.3...(t-3) pr« 1 

2 6 |2 • 16 "6" ■" 2A...i(i+2)'2A..{i-2)i *b' '"]' 

or 
ZT.ttJI^o^Ip^^ l-3...(t-l) 1.3...(t~3) V I 
2^ \2r 16 V ••• 2.4...;(i+2)'2.4...(i-2)i V^^ '•*]' 
according as the attracted point is internal or external to the 
spherical shell, * being any even integer. All these expres- 
sions may be obtained in terms of surface density, by writing, 
instead of pU, 47rcV. 

10. We may next proceed to shew how the potential of 
a spherical shell of finite thickness, whose density is any solid 
zonal harmonic, may be determined. Suppose, for instance, 
that we have a shell of. external radius a, and internal radius 
a', whose density, at the distance c from the centre, is 

Ti Pfi^i h being any line of constant length. 

Dividing the sphere into concentric thin spherical shells, 
of thickness rfo, the potential of any one of these shells, of 
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radius c, at an internal point distant r from the centre will 
be obtained by writing c for 6, ^ for C, 4^Tr(?do for IT, iu 
the first result of Art. 6. This gives 

To obtain the potential of the whole shell, we must inte- 
grate this expression, with respect to c, between the limits 
a and a. This gives 

Again, the potential of the shell of radius o, at an external 
point, will be 

p 4iTrc^ dc PjC^ 47r p_pC^^y 

h* 2i + l 1^ ^"^ 2i + lh'-^'r'^''^^' 

Integrating as before, we obtain for the pptential of the 
whole shell, 

47r p p (a'^°-a^'*^0 

{2i + l){2i+S)h' * r*^' 

Suppose now that we wish to obtain the potential of the 
whole shell at a point forming a part of its mass, distant r 
from the centre. We shall obtain this by consideriug sepa- 
rately the two shells into which it may be divided, the- 
external radius of the one, and the internal radius of the 
other, being each, r. Writing r for a', in the first of the fore- 
going results, we obtain 

2i + l U ^^ ^^ ' 
And writing r for a in the other result, we obtain 
___47r___ p^ r^-a^ 

(2i + l)(2i + 3) hf r**' • 

Adding these, we get for the potential of the whole 
sphere 

47r pP,( a'--r^ , r«^»-a«^M 

2i+l h' I 2 ^"^ {2i + 3)r*^Y 
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It Is hardly necessary to observe that the corresponding 
results for a solid sphere may be obtained from the foregoing, 
by putting a = 0. 

If the density, instead of being ^ JP^c*, be ^Ptc'*^ similar 

reasoning will give us, for the potential of the thin shell of 
radius c and thickness dc at an internal and external point 
respectively, 



*^ A P r^c'^^dc and ^^ P- P - de. 

And, integrating as before, we obtain for the potential of 
the whole shelly 



(2i + l)('m + i + 3)i&^^*" ?^ ^* ^"^ ^^^^"^^^ P^'"'^- 



47r jo^ a^«+» - -'♦»"«+» 

And, at a point forming a part of the mass, 

47r pP, / a"*-*-^-r""*^ , r"^^^^-a "^^^ 1 \ 
2i + l A"* V m-i + 2 ^ ■*■ m + i + 3 r'-'V " 

11. Suppose, for example, that we wish to determine, in 
each of the three cases, the potential of a spherical shell 
whose external and internal radii are a, a', respectively, and 
whose density varies as the square of the distance from a 
diametral plane. 

Taking this plane as that of xy, the density may be ex- 
pressed by ^^^ or pcV*. Now /a^ = — ^ — . Hence the 
density of this sphere may be expressed as 

The several potentials due to the former term will be, 

2 
writing 2 for i and multiplying by ^, 
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15 "F ^* "" ''^' 105 ^-^^ t' ' 15 A" V 2 *^^ 7r» y- 

And for the latter term, writing for i, and 2 for m, and 
multiplying by K, 



4^ 
12 



P u* „'^ ^ P "'-"" *5 P- ( a*-r* r'-a") 
-j^Aa-a)' 15 A" r ' 3 /*' I, 4 ^ br )• 

And, since P^r^ = — ^ — , we get for the potential at an 
internal point 

at an external point 

P f^TT a'-d^ _^ ,r, , 47r a°-a" ) . 

;?lio5-?- (^^-'^+15-^1' 

at a point forming a part of the mass 

A4i5-(-^ + -7i^j(^^ -^ + Tl-4- + -5rJj • 



12. We may now prove that by means of an infinite series 
of zonal harmonics we may express any function of /Lt what- 
ever, even a discontinuous function. Suppose, for instance, 
that we wish to express a function which shall be equal to 
A from /A = l to /a = X, and to B from /Lt=X to /a = — 1. 
Consider what will be the potential of a spherical shell, 
radius c, of uniform thickness, whose density is equal to A 
for the part corresponding to values of fi between 1 and X, 
and to B for the part corresponding to values of jm between X 
and — 1. 

Divide the shell, as before, into indefinitely narrow strips 
by parallel planes, the distance between any two successive 
planes being cdii. 



Digitized by VjOOQ IC 



62 APPLICATION OP ZONAL HARMONICS 

We have then, for the potential of such a sphere at any 
point of the axis, distant z &om the centre, 

for the first part of the sphere 






^irBeio] ^^ 



and for the latter part 

,r 

These are respectively equal to 

?^/:(p..i'.^p.^....p.^....)* 

at an internal point ; and to 



z Jk 






z 

at an external point. 

Now it follows from Chap. IL (Art 23) that if i be any- 
positive integer, 

whence, since / P,d/A = 0, it follows that 
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Also j P,dfi^l-\ j P,dfjL-=l+X. 

Hence the above expressions severally become : 
For the potential at an internal point on the axis 

M«? ]^A (1-X) + £(1 +X) -^{P.(X)-P.(X)} f 

and for the potential at an external point on the axis 

-^{P,(X)-P,(X)}J-... 

Hence the potentials at a point situated anywhere are 
respectively 

3^[{^(l-X) + £(l + X)}P» 

-^{PM-PM^ 
-^{P.i^)-PM^-- 

at an internal point; 
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and 




27rc'Sc {4(l-X) + £(l + X)}^^^^^ 




-^f[PM-PM''f' 




-'';^m)-PM''-f- 


♦• 


-ir^fiPM-PMf'jSf- 


1 
••"J 



at an external point. 

Now, if we inquire what will be the potential for the 
following distribution of density, 

K^(l-X)+-B(l + X)-(^-5){P.(X)-Po(X)lP» 

-(^-P){P.(X)-P.(X)}P.W-... 
-{A-B)[P^,{\) -P,.(X)}P,(^)-...], 

we see by Art. 6 that it will be exactly the same, both at 
an internal and for an external point, as that above in- 
vestigated for the shell made up of two parts, whose densities 
are A and £ respectively. 

But it is known that there is one, and only one, dis- 
tribution of attracting matter over a given surface, which 
will produce a specified potential at every point, both ex- 
ternal and internal. Hence the above expression must 
represent exactly the same distribution of density. That is, 
writing the above series in a slightly different form, the 
expression 

+ {P,(X)-P,(X)}PM + - 
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is equal to A, for all values of/* from J to \ an4 to S for all 
values of fi from \ to — 1. 

13. By a sim^fy^ process, any other discontinuous function, 
whose values ^ve giyen for all values of /i» from 1 to — 1, may 
be expressed. Suppose, for instance, we wish to express a 
function which is equal to A from fi = htQ fjL=^^^, to B from 
/A = X* to /A = Xjj, and to G frorji /* =\, to /i == — ;i. This wiU 
be obtained by adding the two series 

+ {P^.(\)--P^A)}^.W + •••]. 

^^- ^\\+ {-p^w - -p. (\)}^xw + V. 

For the former is equal to -4 — -B from /i = l to fi = \, 
and to from /a = \ to /a = — 1 ; and the latter Is equal to 
JSfrom /i = l to/i = lj, and to (7from /A.=:\,tof6=— 1. 

By supposing A and C each = 0, and 5=1, we deduce a 
series which is equal to 1 for all values of /t from fi = \to 
fi^\ and zero for all other values. This will be 

This may be verified by direct investigation of the 
potential of the portion of a homogeneous spherical shell, 
of density unity, comprised between two parallel planes, 
distant respectively c\ and cK^ from the centre of the 
spherical shell. 

14. In the case in which \ and X, are indefinitely nearly 
equal to each other, letXg = X, andXj=X + dX, We then 
have, ultimately, 

F. H, 5 
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Hence P^.(\0 -Pu,(\) --Pi-X^) - il.(\) 

Hence the series 

+ (2z-M)P,(\)P,(/.) + ...} 

is equal to I'when fp=\ (or, more strietly, when fi has any 
value from X to X -f dX^ and is eqticd to'O for all other values 

of /Lt. 

We hence infer that 

1 + 3P,(\)P»+ ... 4*^(2; + P)P,(\)P//i) + ...- 

is infinite when /i = \, akid z^ro'fcr all other' values of /t. 

15. Representing^ the series • 

. J{l + 3P,(\)P,0^-H,.. + (2i+l)P^)P^tt) + ...} 

by <^(\) for the moment, we see tha* p<fe{\)dK is equal to p 
when /Lt=sX, and to aero for all other values. Hence the 
expression 

is equal to />, when /t* = Xp to p, when /i = \... Supposing 
now that X^, X,... are a series of values varying continuously 
from 1 to — 1, we see that this expression becomes 

j p<l>{\)dK 

p being any function of X, continuous or discontinuous. 
Hence, writing ^(X) at length, we see that 

i|£ />(?X + 8P»|' pP,{\)dK+... 

+ (2z + 1)P,(a.) J^^ /oP,(X)eZX + ...| 

is equal, for all yalues of /i from — 1 to + 1, to the same 
function of fi that p is of X. 
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16. The same conclusion may be arrived at as follows : 

The potential of a spherical shell, whose density is p, 
and volume U, at any point on the axis of z, is 






pdK 



whidiisequalto^jl pd;\+-| /oPj(\)d!\+ ..• 

for an internal point, 

and to 2ftr/°^ + |arP^i(^M^+-^ 

for an external point. 

It hence follows that the potential, at a point situated 
anywhere, is 

for an internal point, 

and to l^ji j" pdK+^^f pP,OC)ctK + ... 

+-;fe^-7V,(x)^+...}. 

for an external point. 

And these expressions are respectively equal to those 
for the potentials, at an internal and external point re- 
spectively, for matter distributed according to the following 
law of density : 

5—2 
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+ (2i + l)P,(ji)j' pP^(K)d\+ .. j. 

It will be observed, in applying this formula, that if p be 
a discontinuous function of X, each of the expressions of the 

form 1 pP/pC)dX will be the jsum of the results of a series of 

integrations, each integration being taken through a series of 
values of X, for which p varies continuously. 
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CHAPTER IV. 

SPHERICAL HARMONICS IN GENERAL. TESSERAL AND SEC- 
TORIAL HARMONICS. ZONAL HARMONICS WITH THEIR 
AXIS IN ANT POSITION. POTENTIAL OF A SOLID NEARLY 
SPHERICAL IN FORM. 

1. We have hitherto discussed those solutions of the 
equation V*F= which are symmetrical about the axis of z, 
or in other words, those solutions of the equivalent equation in 
polar co-ordinates which are independent of ^. We propose, 
in the present Chapter, to consider the forms of spherical 
harmonics in general, understanding by a Solid Spherical 
Harmonic of the i"^ degree a rational integral homogeneous 
function of x, y, Zy of the i^ degree which satisfies the equa- 
tion V'FsO, and by a Surface Spherical Harmonic of the 
i^ degree the quotient obtained by dividing a Solid Sphe- 
rical Harmonic by (a^ +^-f- isO'« Sudi an expression, as we 
see by writing a = rsind cos<^, y = rsindsin^, « = rcosd, 
will be of the i^ degree in sin 6 cos <^, sin sin <^, cos 0\ and 
wiU satisfy the differential equation in Y^ 

or, writing /Lt for cos Oy 

It will be convenient, before proceeding to investigate the 
algebraical forms of these expressions, to discuss some of 
their simpler physical properties. 

2. We will then proceed to shew how spherical har- 
monics may be employed to determine the potential, and 
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consequently the attraction, of a spherical shell of indefinitely 
small thickness* 

We will first establish an important theorem, connecting 
the potential of such a shell on an external point with that 
on a corresponding internal point. The theorem is as follows ; 

If O be tiie centre of such a shell, c its radius, P any in- 
ternal point, F an external point, so situated that P' lies on 
OV produced, and ihatOY .OY ^c\ andifOV^x, OF = r', 
then the potential of the shell at P is to its potential at F 
as c to T, or {which is the same thing) ast' to g. 

For, let A be the point where OP' meets the surface of 
the sphere, Q any other point of its surface. Then^ by a 
known geometrical theorem, 

QP : QP :: AP : AF r.C'-rx^ -d 



And 



r' — o^rr — cr~c* — cr""ci""i^* 



Again, considering the element of the shell in the im- 
mediate neighbourhood of Q, its potential at P is to its 
potential at P' as QP is to QP, that is, as c to r, or (which 
IS the same thing) as r' to c,. which ratio, being independent 
of the position of Q, must be true for every element of the 
spherical shell, and therefore for the whole shell. Hence 
the proposition is proved. 

3. Now, suppose the law of density of the shell to be 

such that its potential at any internal pomt is F(ji,<f>) —^. 

c 

r* 

Then F{fi, ^) --^ must be a solid harmonic of the degree t. 
c 

Hence F(/i, ^) must be a surface harmonic of the degree t. 

Let us represent it by Y^. 

By the proposition just proved, the potential at any 
external point, distant r from the centre, must be 

Y -— 
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Hence, the component of the attraction of the sphere on 
the internal point measured in the direction fropa the point 
inwards, ie. towards the centre of the sphere, is 

And the component in the same direction of the attraction 
on the external point, measured inwards, is 

Now suppose the two points to lie ou the same line 
passing through the centre of the sphere, ajid to be both 
indefinitely close to the surface of the sphere, so that r and / 
are each indefinitely nearly equal to c. 

And the attraction on the external point exceeds the 
attraction on the internal jpoint by 



(2i + l) 



c 



Now, supposing the shell to be divided into two parts, 
by a plane passing through the internal point perpendicular 
to the line joining it with the centre, we see that the at- 
traction of the larger part of the shell on the two points will 
be ultimately the same, while the component attractions of 
the smaller portions, in the direction above considered, will 
be equal in magnitude and opposite in direction. Hence the 

Y 

difference between these components, viz. (2i + l) — , will be 

c 

equal to twice the component attraction of the smaller 

portion in the direction of the line joining the two points. 

But if p^ be the density of the shell, 8c its thickness, this 

component attraction is 27r/>,Sc- 

Y 

Hence (2i + 1) — * = 4t^^p^ Sc, 

c 

2i+l ^ 
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And, if o-^ be the corresponding 8ur£EU^ density, 

It hence follows that if the potential of a spherical shell, 
of indefinitely small thickness, he a surface harmonic, its 
potential at any internal point will he proportional to the 
corresponding solid harmonic of positive degree, and its po- 
tential at any external point will he proportional to the 
corresponding solid harmonic of negative degree. 

That is, the proposition proved for zonal harmonics in 
Chap. IIL Art. 6, is now extended to spherical harmonics in 
general. 

4. The spherical harmonic of th& degi^ee i will involve 
2i + 1 arbitrary constants. 

For the solid spherical harmonic. r*Y], being> a rational 
integral function of x, y, z of the ir^ degree, will consist of 

2 terms. Now the expression V'F, being a 

rational integral function of x, y, z of the degree t — 2, will 

consist of ^^ — ^^ terms ; and the condition that' it must be 

= for all values of x, y, z^ will give rise to'-^ — x-^ relations 

among the ^ — X^ ^ coeflScients of these terms, leaving 

~ ~k — 9"^' ^^ ^** "*" ^^ iiiidependent coeflScients. 

61 We proceed to sheW how the spherical harmtoic of the 
degree i may be arranged in a series of terms, each of which 
majr be deduced by differentiation from the Zonal Harmonic 
symmetrical about the axis of z. The solid zonal harmonic, 
which, in accordance with the notation already employed, is 
represented by r^P^ (/*), iS' a-fanction of i».and r of the degree i, 

satisfying the equation V*F= 0, or -j^ + -7-3- + -7-7 = 0. 
Now, if we denote this Expression by P, {z), we see that 
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since it is a function of z and r, it is a function of the dis- 
tance («) from a certain plane passing through the origin, and 
of the distance (r) from the origin. Further, if we write for z 
the distance from any other plane passing through the origin, 

d^V cPV cPV 
leaving r imaltered, the equation -i-j + XT + "3~8 =0 will 

continue to be satisfied"/ 

Now «Hra(aj + V"-ly)* « being anjr quantity whatever, 
represents tlte distance from a certain plJEine passing through 
the origin^ since in this expression, the suila of the squares 
of the coeflScients of z, x, y is equal to tmity. Hence 
P^[« + a(a? + >/— ly)} is a solid zonal harmonic of the 

degree t, its axis being the imaginary line -- == ^ — = z. 
Therefore the equatioli 

da?^ djt dz""^' 
is satisfied by V^T^ {« + a (a? + V— ly)}, that is, expanding 
by Taylor's Theorem, it is satisfied by 

"*" 1.2.'.. i dz' ' 

for all values of cu 

Hence, since the equation in V is linear, it follows that 
it is satisfied by each term separately, dr that, besides P^ {z) 
itself, each of the * expressions, 

(.+vzij,)4W.(,+vzr,).^,...(»w=Tj,)-*^>, 

satisfies the equation F=isO. 

By similar reasoning we may shew that each of the i ex- 
pressions, 

(»-V^i„*^), (^^i„.^W ...(.-V^,).^, 

satisfies the same equation. 
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Now each of the 2% solutions, thus obtained, is imaginary. 
But the sum of any two or more of them, or the result 
obtained by multiplying any two or more by any arbitrary 
quantities, and adding the results together, will also be a 
solution of the equatioiL Hence, adding each term of the 
first series to the corresponding term of the second, we ob- 
tain a series of i real solutions of the equation. Another 
such series may be obtained by subtracting each term of the 
second series from the corresponding term of the first, and 
dividing by V— 1. We have thus obtained (including the 
original term P<(«)) a series of 2i + l independent solutions 
of the given equation, which will be the 2i + 1 independent 
solid harmonics of the degree %. 

6. We may deduce the surface harmonics from these by 
writing r sin ^ cos ^ for a?, r sin sin (f> for y, r cos for z, 
and dividing by r*. Then, putting coaOssfi, and observing 

that P,{z)^r'P,(ji), ^^ = r*^:^M.... we obtain the fol- 

lowing series of 2i+ 1 solutions : 

cos^sin^^SW, cos2^sm»^^^^, .., cos t> sin*^ ^^ , 

sin^ Sin '^ ' , Sin 2^ sm*^ ~~if* ' • • • ^^ *^ sin'^ *Y • 
Expressions of the form 

(7co8<r^sm'^^^^\ 

or 8Bm<r<l>Bux<'d^^^ , 

or their equivalents, 

58in«r^(l-My^^, 
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(0 and 8 denoting any quantities independent of and (p) 
are called Tesseral Surface Harmonics of the degree i and 
order a. The particular forms assumed by them when 
a = i are called Sectorial Surface Harmonics of the degree t. 

It will be observed that, since * t is a numerical constant, 

Sectorial Harmonics only involve in the form 

sitf0,OT(l-fl*)^. 

The product obtained by multiplying a Tesseral or 
Sectorial Surface Harmonic of the degree i by r* (that is, 
the expression directly obtained in Art. 5) is called a Tesseral 
or Sectorial Solid Harmonic of the degree i. 

7. We shall denote the factor of a Tesseral or Sectorial 
Harmonic which is a function of 0. that is sin'^ — ~-^ , or 

- d^P (a) 
(1 — /ct')" ^ r i ^y *^® symbol T/'>, or, when it is necessary 



IsLri! 
>(/^)or2',(')(cos^ 



to particularize the qoantity of which it is a function, by 



It will be convenient, for the purpose of comparison with 
the forms of Tesseral Harmonics given in the Micanique 
CSleste, and elsewhere, to obtain TJ^'^'> in a completely de- 
veloped form; 

Now, smce Pi\M = ai ^ a o — • — ^^ — - % we see that 
d/j.' "^ 2*. 1.2.3...* rf/t*+' 

- 2\i.2.s. ..id^'T "i''^ "^"172"'* "•••;• 

= 24 (2i- !)...(»- <r+l)^<-» 
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- J (2i - 2) (2t - 3)...(i - o- - 1) /**-'-« 

+ ^TT^ (2t-4) (2t-5)...(i-<r-3y-'-« 



r_o._i)(t--o--2)(t 
2.4.(2i-l)(2t-3) 






And therefore 

^* 2M.2.3...4 ^^ '^^ r 2(2t-l) '^ 

(^-^)(^•-,r-l)(^•-«r-2)(»-«r-3) ,.,_, ), 
■*" 2,4(2i-l)(2t-3) '^ "■/' 

The form given by Laplace for a Tesseral Surface Har- 
monic of the degree i and order a is (see Mecanicpjia Celeste, 
liv. 3, C!hap. 2, pp. 40—47) 

^ (1 -m")' 1^*-- S:i^^^> ^<— « + ...j COB a,^, 

A being a quantity independent of and ^ The factor of 
this, involving fi, is denoted by Thomson and Tait {Natural 
Philosophy, Vol. 1, p. 149) by the symbol 0/*^>. Thomson 
and Tait also employ a symbol ^/•'>, adopted by Maxwell in 
his Treatise on Electricity and Magnetism, Vol. 1, p. 164, 
which is equal to 

(« + o-)(» + <r-l)...(i-<r + l)^ '^/ d/t' * 
or 2' 1.2.-& j,M 
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Heine represents the expression 

(^-,cr)(eW-l)(i-»-2)(t-cr-3) 1 

■^ 2.4.(2*-l)(2i-3) ^ "Y 



or (—1)' ®/*^^ by the symbol PJ(jj), and calls these expres- 
sions by the name Zugeordnete Functionen Erster Art {Hand- 
huch der Kugelfimctionen, pp, 117, 118) which Todhunter 
translates by the term "Associated Functions of the First 
Kind/' which w« shall adopt 

Heine also represents the series 

(t-o-)(;-cr-l) , J 
^ 2(2i-l) ^ 

^ (i-cr)(^-<r-l)(t-<r-2)(t-<r-3) , 

"^ 2.4(2i-l)(2t-3) '^ 

by the symbol ^^(/i), (p. 117). 

The several expressions, Ti^\ @i^\ ^^), Pt, ^i, are con- 
nected together as follows ; 

2M.2.3...i 



;(2i-l)...(i-cr + l) 



r/') = e,( 



.w 



-(i+<r + l)(i + cr + 2)...2r^ ^ ^^ "^^ ^^ ^^ ^^• 

8. It has been already remarked that the roots of the 
equation P, = are all real. It follows also that those of the 

equations ^= 0, -j-i ?= 0.,.^are real also. Hence we may 

arrive at the following conclusions, concerning the curves, 
traced on a sphere, which result from our putting any one 
of these series of spherical harmonics = 0. 

By putting a zonal harmonic=0, we obtain i small circles, 
whose planes are parallel to one another, perpendicular to 
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the axis of the zonal harmonic^ and Bjrmmetrically situated 
•with respect to the diametral plane, perpendicular to this 
axis. If i be an odd number this diametral plane itself 
becomes one of the series. 

By putting the tesseral harmonic of the order cr=0, we 
obtain i — cr small circles, situated as before, and & great 
circles, determined by the equation coscr^ = 0, or sin a^ = 0, 
as the case may be, their planes all intersecting in the axis 
of the system of harmonica^ the angle between the planes of 

• IT 

any two consecutive great circles being - . 

By putting the sectorial harmonic —O, we obtain i 
great circles, whose planes all intersect in the axis of the 
system, the angle between any two consecutive planes being 

IT 

% 

9. The tesseral harmonic may be regarded from another 
point of view. Suppose it is required to determine a solid 
harmonic of the degree t, and of the form F^r*, such that Y^ 
shall be the product of a function of /t, and of a function of ^, 
which functions w^ will denote by the symbols M^, <l>^, respec- 
tively. The differential equation, .to which this will lead, is 

i(.>i,2f,*.+|{(i-.-)f}*,+li.^-.o. 

Now this will be satisfied, if we make Jk^ and <i>^ satisfy 
the following two equations : 

The latter equation gives 

^< = (7 cos (T^ + (7 sin tr^. 

, . Aud, taking o' as an integer, positive or negative, the 
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former is satisfied by M,^ Ty\ i.a (l-,^«)if^j (1 -/*')*, 
as we proceed to prove. 

We know that 

Differentiate o- times, and we get 
whence, by Leibnitz's Theorem, 

or 

^(r+2p /7»'+iP rfi^'P 

and, multiplying by (1 — /tt')^, 

. +(,--<r)(t-+<r+l)(l-.;»f^' = 0...(l). 

Now, putting (l-^'/p-r/'). 
we get 
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-2(.+l)M(l-/*V^'-<r{(l-^')^-<r^«(l-,*V-}p 

- [<r (a + 1) (1 - M'Y- o'il-r /*')^" j0* • 
And i{i+l) ?;(') = »(» + 1)(1-^«)^|5. 

Z^. ^+2 p €. /7<r+l p 



rf/.-"'' ^^ '•^ -di? 



l-M^ 



9 -^1 • 



Hence the equation above given for M^ is satisfied by 
Jfj = T/*'^ and the equation in Y^ is satisfied by 

r,=:(72;<<^>coscr^+ (7'2'/')sincr^. 

10. In Chap. n. Art. 10 we have established the fundamental 
property of Zonal Harmonics, that if i and m be two unequal 

positive integers, I P^PJlfi = 0. This is a particular case 

of the general theorem that if 1^, Y„^ be two surface har- 
monics of the degrees i and m respectively, 



j\j'^T,YJ,^ = Q, 
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For, let Vi, V^ be the corresponding solid harmonics, so 
that V,=r'Y^, V^^r'^Y^. Then, by the fundamental pro- 
perty of potential functions, we have at every point at which 
no attracting matter is situated, 

da?^ dit dz* ' da? ^ df ^ d^ ~^' 
and therefore 

or, in accordance with our notation, Flv*'^«"" Ff»V*^ = 0. 

Now, integrate this expression throughout the whole 
space comprised within a sphere whose centre is the origin 
and radius a, a being so chosen that this sphere contains no 
attracting matter. We then have 

///(F^vT,- F„v'F,) dxdydz^O. 

But also, when the integration extends over all space 
comprised within any closed surface, we have 

///(F.V-F.- V^^JT) dxdydz^jj{v/-^- V„^)d8=0, 

dS denoting an element of the bounding surface, and -r- 

differentiation in the direction of the normal at any point. 

Now, in the present case, the bounding surface being a 
sphere of radius a, and FJ, F^ homogeneous functions of the 
degrees t, m, respectively, 

and, the integration being extended all over the surface of 
the sphere, the limits of /* are — 1 and 1, those of ^, and Stt. 
Hence 



P.H. 
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whence, ifm — ibenot^O, 

The value of / I Y^dfidj^ will be investigated here- 
after. 

11. We may hence prove that if a fwn/cUon of ii and ^ 
can be developed in a series of surface harmonics, such de- 
velopment is possible in only one way. 

For suppose, if possiblei that there are two such develo{>- 
ments, so that 

and also 

Then subtracting, we have 
0= r^- F;+ Fj- F;+ ... J^T.^Ti-^- ... identlcaUy. 

Now, each of the expressions F, — F/, F^ — F/... F^ — Yl 
being the difference of two surface harmonics of the degree 
0, 1, ...t... is itself a surfetce harmonic of the degree 
0, 1, ...t.... Denote these expressions for shortness by 
^, ^j ... ^ ... so that 

=4+ -^^+ ...+ ^< + ... identically. 

Then, multiplying bv Z^ and integrating all over the 
surfiEtce of the sphere, we have 



0=r r Z^dfid^ 



That is, the sum of an infinite number of essentially 
positive quantities is =0. This can only take place when 
each of the quantities is separately = 0. Hence Z^ is identi- 
cally = 0, or Yl «= Y^ and therefore the two developments 
are identical 

We have not assumed here that such a development is 
always possible. That it is so, will be shewn hereafter. 



Digitized by VjOOQ IC 



TE8SEBAL AND SECTORIAL HABMONICS, 83 

12. By referring to the expression for a surface har- 
monic given in Art 4, We see that each of the Tesseral and 
Sectorial Harmonics involves (1 — ft')*, or some power of 
(1 — ij!)\ as a factor, and therefore is equal to when /tt = ± 1. 
From this it follows that when /a = ± 1, the value of the 
Surface Harmonic is independent of ^, or that i£T{ji,^) repre- 
sent a general surface harmonic, Y [± 1, <^) is independent of 
ff>y and may therefore be written as F (± 1). Or F (1) is the 
value of x{ji, <f>) at the pole of the zonal harmonic Pii/J^), 
r'(— 1) at the other extremity of the axis of JP^ (/a). 

We may now prove that 



/^ 



r^=2,rr,(i)P,o*). 



For, recurring to the fundamental equation, 

Now, if we integrate this equation with respect to ^, 
between the limits and 27r, we see that, since 



/ 



d^'Y dY. 



and the value of Y^ only involves <f> under the form of cosines 
or sines of ^ and its multiples, and therefore the values of 

-^' are the same at both limits, it follows that 



L 






Hence 



|i{(i-/*')(/rw)}+*(i+i)(f;w)=o. 

Hence I Y^d(|> is a function of fi which satisfies the 

fundamental equation for a zonal harmonic^ and we therefore 
have 



6—2 
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r^ 



G being a constant, as yet unknown. 

To determine O, put /*=!, then by the remark just made, 
r^ becomes F,(I), and is independent of ^. Hence, when 

/t=l, pr,(?^-2Trr,(l). AlsoP,(/t) = l. We have there- 
fore * 2irr,(l) = 0; 

It follows from this that 

13. We may now enquire what will be the value of 

j[J^Y,Z,d,jLAf>, 

Y^, Z^ being two general surface harmonics of the degree i. 
Suppose each to be arranged in a series consisting of the 
zonal harmonic P^ whose axis is the axis of z, and the system 
of tesseral and sectorial harmonics deduced from it. Let us 
represent them as follows : 

r,= AP, 

+ G^Tp^^ cos ^ + C^T^^^ cos 2^ + ... + C^T/') cos o-^ + ... 

+ Cir«cosi^ 
+ )8;r/«sin<^ + i&,r/2)sin2^+... + i8;r/*^>sincr<^ + ... 

+ i8ii;Wsint^; 

+ c^T^Ci) cos ^ + oj:'l^'^ cos 2^ + ... + c,r/'> cos o-^ + . .. 
+5j3;ci)sin^ + «j^r/^ sin 2^ + ... +5^7/*^) sino-^ + ... 

Hence the product Y^Z^ will consist of a series of terms, 
in which j> will enter under the form cos <r^ cos a-'<^, or 
cos cr^ sin o"'^. This expression when integrated between 
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the limits and 27r vanishes in all cases, except when 
a^a- and the expression consequently becomes equal to 
cos'o-<^, or sin\o-^. In these cases we know that, a being any 
positive iAteger, 



/•2ir TBir 

I co^a^d(f>^\ siv? a(l>d(}) = 
Jo J-o 



TT. 



Hence the question is reduced to the determination of the 
value of 

-1 



/: 



Now i;w=(i-m')^^ 

2'. 1.2. 3.,.*^ '^^ d(>.^^' * 

But, by the theorem of Rodrigues, proved in Chap. II. 
Art 8, we know that 

Hence T/') may also be expressed under tiie ^orm 

^ ^^ 2M.2.3:..iL^^*'^ ^^ d/i*-' ' 

whence it follows that 

^^« i -^ "-J V2',1.2.3...i7 \i-a - dfjf+' d/i^-' ' 

Now, putting (ji' — iy^M for the moment, and inte- 
grating by parts, 



/ 






d/A»+' rf/4*-' '^ dfJl.*'^''-^ rf/4*-' 



-/ 
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The factor -j—. — vanishes at both limits, hence 

by a repetition of the same process. 

And by repeating this process o- times, we see that 

= (- 1)' (2«. 1.2,3 ...»)'|* P*diJ. 
= (-l)'(2'.1.2.3...»)'^. 

Hence /^W<?m = ^2?T1' 

and therefore 

r r{T,^'^ COS <r(l>ydfid<f>r.p ^ {TifHmafi'dfid(t> 

It will be observed that this result does not hold when 
cr =x 0, in which case we have 

Hence f P^YtZ^diid^ 

Afl, 



2i + l 



* In this case f cos* <r<l>d<l>-j sin' (r0(2^=:2r. 
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+fef (^<^<'+^^')+•••+L?!(^A+^i«l)^ 



l4 We have hitherto considered the Zonal Harmonic 
under its simplest form, that of a " Legendre's Coefficient ** in 
which the axis of «, ie. the line from which is measured, is 
the axis of the system. We shall now proceed to consider it 
imder the more general form of a "Laplace's Coefficient,'* 
in which the axis of the system of zonal harmonics is in any 
position whatever, and shall shew how this general form may 
be expressed in terms of PiQi) and of the system of Tesseral 
and Sectorial Harmonics deduced from it. 

Suppose that ^, (f/ are the angular co-ordinates of the 
axis of the Zonal Harmonic, i.e. that the angle between this 
axis and the axis of z is 6^, and that the plane containing 
these two axes is indined to a fixed plane through the axis 
of z which we may consider as that of zx, at the angle ^\ 
In accordance with the notation already employed, we shall 
represent cos ^ by fjf. 

The rectangular equations of the axis of this system 
will be 

.« y 



sin ff cos ^' sin ff sm ^' cos ff ' 

Hence the Solid Zonal Harmonic of which this is the axis 
is deduced from the ordinary form of the solid zonal har- 
monic expressed as a function of z and r by writing, in place 
otz, xsinffcoB<f>' + ymiffsm<f>+zcosff. 

To deduce the Surface Zonal Harmonic, transform the solid 
zonal harmonic to polar co-ordinates, by writing r sin 5 cos 
for a?, r sin 5 sin ^ for y, r cos for z, and divide by r*. 

The transformation from the special to the general 
form of surface zonal harmonic may be at once effected, 
by substituting for fi, or cos 0, cos^cosa +sin^sin^cos(^— ^'). 

Now, in order to develope 

Pi {cos ^cos ^ + sin ^ sin ^ cos (^ — ^')} 
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in the maimer already pointed out, assume 

P, {cos Ocosff + sin sin ff cos {(f> — (j>)} 

' =^P^0A) + ((7Wcos<^ + /SfWsin^) 3^(1) 

+ (0(2) cos2<^+ /8f(2) sin 2^) i;(2) + ... 

+ (0<'>cos 0-^+ /8f<<^>sin a^ T^^^H... 

+ {G(^ cos i<l> + 8^hmi<l>) Ti% 

the letters -4, ... (7^*^^, 8^^K.. denoting functions of /*' and 
^', to be determined. 

To determine C^^\ multiply both sides of this equation 
by cos a^T}^"^ and integrate all over the surface of the sphere, 
ie. between the limits — 1 and 1 of /Lt, and and Stt of ^.. 
We then get 

r j 'P, {cos e cos ff + sin 5 sin ^ cos (</> - f )} cos cr^2;<<^> diid<f> 
= a<<')r [^ (cos <T<i>T^^)^ d,jd<l> 

""[i-^2t + l^ • 

It remains to find the value of the left-hand member of 
this equation. 

Now cos cr^iy^) is a surface harmonic of the degree i, and 
therefore a function of the kind denoted by Yj in Art. 12. 

And we have shewn, in that Article, that 

that is, that if any surface harmonic of the degree i be multi- 
plied by the zonal harmonic of the same degree, and the product 
integrated all over the surface of the sphere, the integral is 

equal to ^' — T ^'^^ ^^ value which the surface harmonic 

assumes at the pole of the zonal harmonic. 
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Hence 

and therefore 

I I P, {cgs 5 cos ^ + sin 6 ^mff cos (^ — <!>)] cos tr^T\''^dfidj> 

or C<') = 2 [t^cos «rA' I?') OtO- 
I — <r 



Similarly /SfW = 2 h= sin <rf 2;W Oi')- 

And to detenoine A, we have 
Y^ j P, {cos ^ cos ^ + sin ^ sin ^ cos (^ - ^') } P, O) <?/*# 

•'• 2?+l-*«^^~^2?+I' 
or 4 = P>')- 
Hence, P, {cos ^ cos ^ + sin ^ sin ^ cos (^ - <}>')} 

= P.0*')-P.(M)+2^co8(^-f)T,W0*')2;«(/*) 

+ 2L^cos2(^-f)r/»0*')2;<«>M + -.. 
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+ 2 [^cos«r (^-f ) r/') (/tO 2?'> M + - 

15. We have already seen (Chap. ll. Art. 20) how any 
rational integral function of fi can be expressed by a finite 
series of zonal harmonics. We shall now diew how any 
rational integral function of costf, sin tfcos^/sintf sin^, 
can be expressed by a finite series of zonal, tesseral, and 
sectorial harmonics. 

For any power of cos ^ or sin ^, or any product of such 
powers, may be expressed as the sum of a series of terms of 
the form cos a-<f>, or sin a^, the greatest value of <r being the 
sum of the indices of cos <f> and sin <f>, and the other values 
diminishing by 2 in each successive term. Hence any 
rational integral function of cos tf, sin tf cos 0, sin 5 sin <f), will 
consist of a series of terms of the form 

cos** sin* 0coa (r<f> or cos'^ sin* ^sin <r<f), 

where n is not less than <r. 

If n be greater than <Tyn-'<r must be an even integer. Let 
71 — <7 = 25, then writing sin"d under the form (1 — cos* tf)' sin*'^, 
we reduce cos** sin* cos atf> to the sum of a series of terms 
of the form cos' sin*' cos o-^, or, writing cos d = /i, of the 

form fiF (1 — /A*) * cos a^. 

Similarly cos**^sin*dsin<7^ is reduced to a series of 

terms of the form /if (1 — /i-')*sin<r^ 

^''^ ^° (p.,.«.)(p+<.li)...(p+i) ^>^^'> 

and fiP'^*' can be developed in a series of terms of the form 
of multiples of Pp+<r> Pp+<r-2 .... (Chap. n. Art. 17.) 

Hence /^ can be expressed in a series of the form 

^r^ (-^0 Pp-^rv + -^a -^+0^8 +•••)> 
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Aq, A^ representing known numerical constants, and therefore 
fiP {1 -- /Ji^y assumes the form 

consequently multiplying these series by cos<r^ or sin <7^, we 
obtain the developments of 

/aP (1 — /A*) * COS <7^ and /t**^ (1 — /u.*)* sin <7^ 
in series of tesseral harmonics. 

16. We will give two illustrations of this transformation. 

First, suppose it is required to express cos* sin'5 sin ^ cos^ 
in a series of Spherical Harmonics. 

Here we have sin ^ cos ^ = ^ sin 2^. 

Hence cos*5sin*5sin^cos^ = ^cos'dsin*dsin 2^. 

Comparing this with cos** sin* sin o-^, we see that n is 
not greater than <r. 

Hence cos*dsin*^sin^cos^ = ^/x'(l — /A')sin2^. 

A J 1 1 t?" 4 

and ;t« = ^P^ + |p^ + lp^, 



. , 1 fScPP, 4<fPA 
'• ^ ""12135 dfi,'^7 dp.*) 
2 (fP. 1 <f P, 
"l05d/t''*"21d/*»' 
.'. cos* sin* d sin ^ cos ^ 
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Next, let it be required to tranflform cos'd sin' sin <p cos*<^ 
into a series of Spherical Harmonics. 

Here sin ^cos^ ^ = ^ sin 2^ cos ^ = t (sin 3^+ sin ^). 

Now cos' d sin* d sin 3^ = /a' (1 - /a')' sin 3<^ 

Also cosM sin" ^ sin ^ =/*» (1 - /t") (1 -/*»)* sin ^ 
= (/*»-/)(!-/*')* sin ^ 



"Sl"'-?!"')''-"^'-*- 



Also (Chap. n. Art. 17) 

„«_lip^24 10 1 
'^ ~231 ' 77 * 21 * ? •* 
Hence cos' d sin' 6 sin 3^ 

1/16 <PP, , 24(?PA „ i^. . _ , 
= 120(23-1^ +77^)^^-'*^'^"^^'^ 

= {3^5 ^."^ + 3-15 ^«^3^ 

A J »a • ia • J. f 8 dP. , 4! dP. , 5 dP, 
Andcos'5sm'dsm*=-^gg3^« + ^^ + ^^' 



2 dP, 1 dPA ., ^i • . 

/ 8 dP, 2 (iP, 4 dPA ,, ^1 . . 
= -(693d;f-385-^-63^j(^-'*^*^^ 

= -/'— T«- — TW-— T'^Asind' 
W • 385 * 63 « j^^y' 

cos'dsin'flsin^cos«^=|g^2;(») + j^2;(«}sin3^ 

-J— ra)__Ly(i) — 5L r (1)1 sin«*. 
1693 • 770 * 63 * ) ^ 
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17. The process above investigated is probably tbe most 
convenient one when the object is to transform any finite 
algebraical function of cos d, sin d cos 0, and sin sin ^, into 
a^ seizes of spherical harmonics. For general forms of a 
function of ^ and ^, however, this method is inapplicable, 
and we proceed to investigate a process which will apply 
universally, even if the function to be transformed be discon- 
tinuous. 

We must first discuss the following problem. 

To determine the potential of a spherical shell whose 
surface density is jP(/Lt, ^), JP denoting any function whatever 
of finite magnitude, at an external or internal point. 

Let c be the radius of the sphere, r' the distance of the 
point from its centre, ff^ ^' its angular co-ordinates, V the 
potential. Then /a being equal to cos 5 






F{fi,<f>)<?dfid(l) 



[r'^-2cr'{cos5 cose'+ sin5 sin 5' cos (^-^Ol + ^* *- 

The denominator, when expanded in a series of general 
aonal harmonics, or Laplace's coefficients, becomes 

for an internal and an external point respectively, P^ (/a, ^) ' 
being written for 

Pi {cos ^ cos ^ + sin d sin ff cos (^ — ^') }• 

Hence, V^ denoting the potential at an internal, F,, at an 
external, point, 
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+ ... + p?|J*'P, (Ml ^) -Z^0^> ^) ^/^# + — f • 

It will be observed that the expression Pi{fi, <f>) involves 
fi and fi symmetrically, and also <f> and <)>. Hence it satisfies 
the equation 

And^ since fi and 6 are independent of fi and ^\ this 
diflferential equation will continue to be satisfied after P, has 
been multiplied by any function of [a and <)>, and integrated 
with respect to /ju and 0. That is, every expression of the 
form 

is a Spherical Surface Harmonic, or "Laplace's Function" 
with respect to fju and <}) of the degree i. And the several 
terms of the developments of T\ are solid harmonics of the 
degree 0, 1, 2...i... while those of F, are the corresponding 
functions of the degrees —1, —2, — 3... — (i+ 1), ... And 
these are the expressions for the potential at a point (r', fi, <f)) 
of th^ distribution of density F(ja, <f>') at a point (c, fi\ <(>), 

Now, the expressions for the potentials, both external 
and internal, given in the last Article, are precisely the same 
as those for the distribution of matter whose surface density is 

+ {2i + 1) fJ^PiiP'> i>) F(fi, <l>) dfid<f> + ...| , 
or, as it may now be better expressed. 
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+3 r p2?,{co8dco8d'+Bm^sm^oos(<^-<^')^0^i^)d/*^ 






(2i+l) j j P,{cos^cos^+sin5sin^cos(^^')}^(/*>^)^/^^^+-- 

And, since there is only one distribution of density which 
will produce a given potential at every point both external 
and internal, it follows that this series must be identical 
Yfith. F{fj/, <f>y We have thus, therefore, investigated the 
development of F(/i, 0') in a series of spheric^ surface 
harmonics^. 

The only limitation on the generality of the function 
.F(jJby if)) is that it should not become infinite for any pair of 
values comprised between the limits —1 and 1 of /i^ and 
and 27r of ^. 

18. Ex. To express cos 2<^' in a series of spherical har- 
monics. 

For this purpose, it is necessary to determine the value of 
{2i + l)j j 'PJcos^cog^+sin^sin^cos(^-^')}cos2^d^. 

Now Pi {cos ^cos ^ + sin ^ sin d' cos (^ — (f>)] 

= P,(cos^)P,(cos^ 

. .2 . ^dP,(cos0) . ^dP, (cos ff) ,^ .,. 

2 
■^(i-l)t(i + l)(t + 2)®"''^ 

c?P,(cos^) . ,^<?P,(oos^) o.. .^ . 
%^ ^ "^^ %^ ^ cos2(<^^^0 + >» 

Now / oos<r(^-^')oos2^{i^ = 0, 

for all values of <r except 2, 

* In conneotion with the Bubjeot of this Article, see a paper by Mr G. H. 
Darwin in the Mesienger of Mathematics for Maroh, 1877* 
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/•Sir 

And I cos 2 (^ — ^') COB 2^ (2^ SETT cos 2^'. 

JO 

Also 
And 

f <?">*-!)% _, <y(A*'-iy <^"o**-iv 

= 2M.2.3...t/.P.- ^'^^.7^)' . 
Now when /tt == 1, 

And when ft = — 1, 

Hence 

£^Bin'ggc?^ = ^,^j3 . g.l.2.3...z{l-(-in 

= 4 or 0, as i is even or odd ; 

.-. r J" 8in*fi ^^^^cos 2 (0 - f ) cos 2^ dt^d^ 

= 47r cos 2^' or 0, as i is even or odd ; 
.'. cos 2^' 

If- 2 . . ,„d'P,(cos^) 

,_ « 2 , . .-(ppicos^) _,, 
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+ . 



Hence the potential of a spherical shell, of radius c and 
surface density cos 2<f>, will be 

and 

at an internal and external poiQt respectively. 

19. We will now explain the application of Spherical 
Harmonics to the determination of the potential of a homo- 
geneous solid, nearly spherical in form. The following 
investigation is taken from the M^caniqvs Celeste, Liv. iii. 
Chap. II. 

Let r be the radius vector of such a solid, and let 
r = a + a (a^r^H- ttgF, H- ... + a^r, + ...)> 
a being a small quantity, whose square and higher powers 
may be neglected, o^, a^,.,.a^,.. lines of arbitrary length, and 
Y^y Fj,... F,... surface harmonics of the order 1, 2,...V... re- 
spectively. 

The volume of the solid will be ^ Tra', 

For it is equal to 

JoJ-iJq 

= X ira!^, since I / y^ djidcj) = 0, 

for all values of i. 

F. H. 7 
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Again, if the centre of gravity of the solid be taken as 
origin, a^ = 0. 

For if z be the distance of the centre of gravity from the 
plane of xy^ 

Similarly 
- Tra'u; = 4a' a . a^ I I (1 — /it')* cos ^ F^ c?/irrf(^, 

-7ra'y=4a'a. ttj I I (1— /iA')4sin^r'je?/A(?<^. 

Now Yj^ is an expression of the form 

J/a + 5(1-/a'j*cos(^+ C(l-/A»)isin<^, 

and therefore all the expressions x, y, i cannot be equal to 0, 
unless a^ = 0. 

We may therefore, taking the centre of gravity as origin, 
write 

r ^ a + a{a^Y^-\- .,. •¥ a^Y^•\' ...), 

as the equation of the bounding surface of the solid. 

Now this solid may be considered as made up of a homo- 
geneous sphere, radius a, and of a shell, whose thickness is 

a(a,F,+ ... + a,r,+ ...). 

The potential of this shell, at least at points whose least 
distance from it is considerable compared with its thickness, 
will be the same as that of a shell whose thickness is aa, and 
density 



..gn^...^|'r.+,...). 
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p^ being the density of the solid. Therefore the potential, 
for any external point, distant JR from the centre, will be 

The potential at any internal point, distant S from the 
centre, will be made up of the two portions 

I 'Trp.B' + 27rp. (a' - S^ or 2wp (a' - ^^ 

for the homogeneous sphere, 

for the shell, and will therefore be equal to 

20. If the solid, instead of being homogeneous, be made 
up of strata of different densities, the strata being concentric, 
and similar to the bounding surface of the solid, we may 

deduce an expression for its potential as follows. Let :>^ r be 

the radius vector of any stratum, p its density, r having the 
same value as in the last Article, and p being a function 
of c only. Then, he being the mean thickness of the stratum, 
that is the difference between the values of c for its inner 
and outer surfaces, the potential of the stratum at an ex- 
ternal point will be 

'^^i+iBf'^-') ^^^• 

To obtain the potential of the whole solid at an external 
point we must inteofrate this expression with respect to c, 
between the limits and a, remembering that p is a func- 
tion of c. 

7—2 
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Again, the potential of the stratum, above considered, 
at an internal point will be 



47rpcSa + 47rpa — (^-Y»-^+-»7«-^ + ..- 



+ ^ 



aX R' 



j+-') (2). 



2i+l c* 

To obtain the potential of the whole solid at an internal 
point we must integrate the expression (1) with respect to c 
between the limits and R, and the expression (2) with 
respect to c between the limits It and a, remembering in 
both cases that /9 is a function of c, and add the results 
together. 
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CHAPTER V. 

SPHERICAL HARMONICS OF THE SECOND KIND. 

1. We have already seen (Chap. ii. Art* 2) that the 
differential equation of which P, is one solution, being of 
the second order, admits of another solution, viz. 

CPJ ^^ 



^JP,»(1 -/.'')• 



Now if fi between the limits of integration be equal 
to + 1, or to any roots of the equation P^ = (ill of which 
roots lie between 1 and —1), the expression ' under the 
integral sign becomes infinite between the limits .of inte- 
gration. We can therefore only assign an intelligible 
meaning to this integral, by supposing fi to be always be- 
tween 1 and 00 , or between — 1 and — oo • We will adopt 
the former supposition, and if we then put C7= — 1, the 

expression p, ,^ ^ ^v (i.e. p2 / 2 _ -i n ) ^iU he always posi- 
tive. We may therefore define the expression 

^j, p.'(^«-i)' 

as the zonal harmonic of the second kind, which we shall 
denote by ^„ or Q^ (fju), when it is necessary to specify the 
variables of which it is a function. 

It will be observed that, if ^ be greater than 1, P, is 
always positive. Hence, on the same supposition, Q^ is 
always positive. 

We see that 



^ r dfi 1, /A + 1 
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An4 in a similar manner, the values of Q^, Q^,... may 
be calculated. 

2. But there is another manner of arriving at these 
functions, which will enable us to express them, when the 
variable is greater than unity, in a converging series, with- 
out the necessity of integration. 

This we shall do in the following manner. 

Let ?7= , V being not less, and /jl not greater, than 

unity. 

Th ^-. 1 ^ 1 

n^.x,dU_ i;'-l ,. ^dU_ l^u? 
^^ "^^ di^'Jir:^' ^^''f^^ dfjL'^^^^^'' 

Now, let be expanded in a series of zonal harmonics 

P», P»..,P,(/t),sothat 

^=r^ = *,(") -P.W +^» -P» + .-. + ^.M ■?.(/»)+... 

Then|^|(l-/t')^^=...-«(,-+l)^.(v)P.(/*) + ... 
by the definition of P (/x). 
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^ '^ i {('-^ f }=-4 {(I-') T} ^.w *■■■ 

And these two expressions are equal Hence, equating 
the coefficients of P^ (/*), 






i{('-*)^'}=--('-^^)*.w- 



Hence 0,(i/) satisfies the same differential equation as P, 
and Q^. But since U=^ when i/ = oo , it follows that ^, {v) =0 
when !/= 00 . Hence ^^{v) is some multiple of Qi{y)==AQ^{v) 
suppose. It remains to determine A. 

Now, <f)^{v) may be developed in a series proceeding by 
ascending powers of - , as follows. 

We have = -+^ + --+-^+ •••• 

v — fi V ir V 

and also = <^oW^oW +^iW J^iW +.-. + ^iM ^< W+.- 

Now, by Chap. ii. Art. 17, we see that, if m be any 
integer greater than % the coefficient of P^ in /a** is 

/«• . ^N (w — i+ 2) (m-i + 4)...(m-l) .« . , ii 
(^^•^^) (m+\- + l)(m + t-l)...(m4-4)(J+2) ^^^ ''^^' 

and (2t + l) , jTiTf ""5^".; 'T^kr xi. ^^'^^^^^"' 
^ ' (m+»+l) (m + t— l)...(»i+3) (wt+1) 

m — t being always even. 

Hence, writing for m successively i, i + 2, » + 4, ... we get 

./x ,«. ,N f 2.4...(i-l) 1 

^.(0 = (2» + l)|^2i + l)(2i-l)...(i + 2).- 

4.6...Ct + l) 1. 

■*" (2t + 3) (2i + !)...(»• + 4) v**» 

6.8...(i + 3) 1 , \ -t .■> AA 
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2.4...t 1 



= (2»+i){(27^ 



and v-.n^ ^j v/..^ij ^2i-l)...(i4-l) v^'^ 



+ 



4. 6... (1 + 2) 1 



(2f+3j(2» + l)...(i + 3) v' 



,*+» 



6.8...(t + 4) 1 , , .. „ 

+ (2i.H5)(2» + 3)...(» + 5) /^»+ --^ ^f » ^« «^«°- 



..I if t be ( 



Now, recurring to the equation 

•'"■PiW; (a* -1) 



we see that, if ^,(1*) be developed in a series of ascending 

powers of - , the first term will be ^ ... — , , ,^, , where C 

is the coefficient of /** in the development of P, (/*) ; 

,, ^. „ (t+-2)(t + 4)...(2t-l) .. ., ,, 
that IS C=^ — 2.4.6...(i-l) — ^^^^^eodd, 

and = (^->-l)(^^H- 3)^:^4^ )■■■(2^-l) ^^ . ,^ ^^^„ 

2.4. 6...t 

Hence the first term in the development of Q, (v) is 
2.4.6...(t-l) 



(i + 2)(i + 4)...(2i-lj(2i+l) 



if i be odd, 



, 2.4.6...i .. ., 

(i+ 1) (i + 3)...(2i- 1) (2^ + 1) 

which is the same as the first term of the development of 
P, (p), divided by gfipy • 

Hence A — 2i-\- 1, and we have 
1 



I/ — /A 



= <?oMPoW+3(?,MP,(/.) + 5(2,(i;)P,(/.) + ... 



3. The expression for Q, may be thrown into a more 
convenient form, by introducing into the numerator and de- 
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nominator of the coefficient of each term/ the factor neces- 
sary to make the numerator the product of i consecutive 
integers. We shall thus make the denominator the product 
of i consecutive odd integers, and may write 

^ 1.2.3...i 1 3.4.5...(i + 2) 1 

5.6. 7.. .(1 + 4) 1 
"^5.7.9...(2t+5)i;'^"^- 

(2k+l){2k+2)...{i-¥2k) 1 
"^ (2&+ 1) (2& H-3)...(2iH- 2k + l) v'^'' "*" "' 
whether i be odd or even. 

4. We shall not enter into a full discussion of the pro- 
perties of Zonal Harmonics of the Second Kind. They will be 
found very completely treated by Heine, in his Handbuch der 
Kugelfunctionen, We will however, as an example, investi- 
do,,, 
du 



gate the expression for -^ in terms of Q^^^, Q^^^, 



Recurring to the equation 
1 



+ (2i + l)5Q.MP.(/*) + ... 
we see that 

+ (2,-+l)5«,W^ + (2.-+3)«„W^^+.... 
Now we have seen (Chap. il. Art. 22) that 
^1^ = (2/- ») P...(/*) + (2» - 5) P,.(/*) + ... 

Hence ^^^ = (2t + 1) P, (/*) + (2t - 3) P^ (/.) + . .. 
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^^ = (2^ + 9) P^if.) + (2i+ 5) P^W 

+ (2i+l)P,(^) + ... 

And therefore the coefficient of P,(/ti) in the expansion 
- rf 1 . 

of -j IS 

dfi v — fi 
(2t+l) {(2i+3) Q^.(i')+(2i+7) <?.„W + (2i+ll) Q^{v) + ...}. 
Again, 



dvv-ii~ dv ^'^J^" dv 



+ (2i + l)^P.W. 



And -^ +:7 = 0. 

dv v — fi diiv — fj* 

» 

Hence, comparing coefficients of P, (/a), 

-(2;+ii)Q,„w-... 



^^) = -(2i + 3)Q.„W-(2i+7) <?^W 



Hence it follows that 
and therefore that 

5. By similar reasoning to that by which the existence of 
Tesseral Harmonics was established, we may prove that there 
is a system of functions, which may be called Tesseral Har- 
monics of the Second Kind, derived from T^^ in the same 
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manner as Q, is derived from P,. The general type of such 
expressions will be 



^■'"w/. V«X-i) - 



and this when multiplied by cos<r<^or sino-^, will give an 
expression satisfying the diflferential equation 

and which may be called the Tesseral Harmonic of the 
second kind, of the degree i and order c. 
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CHAPTER VI. 

ELUPSOIDAL AND SPHEROIDAL HARMONICS. 



1. The characteristic property of Spherical Harmonics 
is thus stated by Thomson and Tait (p. 400, Art. 537). 

"A spherical harmonic distribution of density on a spheri- 
cal surface produces a similar and similarly placed spherical 
harmonic distribution of potential over every concentric 
spherical surface through space, external and internal." 

The object of the present chapter is to establish the ex- 
istence of certain functions which possess an analogous pro- 
perty for an ellipsoid. They have been treated of by Lamd, 
in his Lefons sur les fonctions inverses des transcendantes et 
les fonctions isothermes, and were virtually introduced by 
Green, in his memoir On the Determination of the Exterior 
and Interior Attra^ctions of Ellipsoids of Variable Densities, 
(Transactions of the Cambridge Philosophical Society, 1835). 
We shall consider them both as functions of the elliptic co- 
ordinates (as Lam^ has done) and also as functions of the 
ordinary rectangular co-ordinates; and after investigating 
some of their more important general properties, shall pro- 
ceed to a more detailed discussion of the forms which they 
assume, when the ellipsoid is a surface of revolution. 

2. For this purpose, it will be necessary to transform 
the equation 

into its equivalent, when the elliptic co-ordinates e, v, v* are 
taken as independent variables. If a, b, c be the semiaxes 
of the ellipsoid, the two sets of independent variables are 
connected by the relations 
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a? y' . g' _, 



Thus a* + e, 5" + e, c' + e are the squares on the semiaxes 
of the confocal ellipsoid passing through the point a?, y, z, 

c? + 1;, 6* + 1;, c' + u, the squares on the semiaxes of the 
confocal hyperboloid of one sheet. 

a' + v\ J' + V, c' + v\ the squares on the semiaxes of the 
confocal hyperboloid of two sheets. 

Thus, € is positive if the point a?, y, z be external to the 
given ellipsoid, negative if it be internal. 

And, if a? be the greatest, c' the least, of the quantities 
a\ V, &, 

€ willlie between — c' and 00 , 



V 


a 


99 


-b* 


99 


-c\ 


V 


» 


99 


-a* 


J> 


-b\ 



3. Now -Tj + -j-i + -j-i = IS the condition that 

taken throughout a certain region of space, should be a mini- 
mum. In the memoir by Green, above referred to, this 
expression is transformed into its equivalent in terms of a 
new system of independent variables, and the methods of the 
Calculus of Variations are then applied to make the resulting 
expression a minimum. We shall adopt a direct mode of 
transformation, as follows : 

Suppose a, fi, 7 to be three functions of x, y, z, such that 

v'a = o, v'/9=0, vV=0 (1), 

such also that the three families of surfaces represented by 
the equations a = constant, /9 = constant, 7 = constant, inter- 
sect each other everywhere at right angles, i.e. such that 
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dxdx dy dy dzdz ' dxdx dydy dzdz" * 

d2dfi^d2d§^d2dff^^ 

dxdx dydy dzdz ^ 

Then 

dV^dVdadVdfidV&y 
dx d% die dfi dx d/y dx^ 

da^^mXdx) ■*" dff'Kdx) "*■ dy" \dx) 

^ ^d'Vdfidy ^ ^d'Vdydx ^ ^d'V dad/3 
d^dydxdx drfdzdxdx dad^dxdx 

^dV^.dVd^ dV^ 
dada^ d^da? dy dm?* 

dW d^V 

-j-2 and -T-g being similarly formed, we see that, when the 

three expressions are added together, the terms involving 

dV dV dV 

-T- , -17T, -7- will disappear by the conditions (1), and those 

i-^volving ^,^^,£^ by the conditions (2). Hence 



^'^ 



Now, let 



=/. 



{(a' + t)(6' + t)(cM-'f)]** 
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^ Ji/ Ua' 



All these expressions satisfy the conditions (1), for a is 
the potential of a homogeneous ellipsoidal shell, of proper 
density, at an external point, and fi and 7 possess the same 
analytical properties. 

Again, a is independent of v and v, and is therefore con- 
stant when € is constant. Similarly /3 is constant when v is 
constant, and 7 is constant when v is constant. Hence a, /8, 
7 satisfy the conditions (2). 

Now 

And a. +j;/r"+-a- 



a* + € V-\-€ C* + € 



i 






with similar expressions for -^ and t- . Hence, squaring 
and adding, 

But from the equations 

a^+e^fc^ + e^c^ + c ' a'+v^ b' -^ v^ c' + v ' 
* ]8 is a purely imaginary quantity. We may, if we please, write J -lp> 
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7 +ijV:7+ 72-7-77 = 1. 



we deduce 

a> being any quantity whatever. For this expression is of 
dimensions in o), €, v, v\ it vanishes when o) = 6, v, or v, 
and for those values of o) only, it becomes infinite when 
6> = — a', — y, or — c*, and for those values of <o only, and it is 
=s 1 when ft) = 00 , 

From this, multiplying by a* + ft), and then putting 
ft) = — a', we deduce 

,_ (€ + a') (i/H-a'Xi/' + a*) 

a result which will be useful hereafter. 

Again, diflferentiating with respect to ft), and then putting 

ft)=:€« 

X* y I '^^ ^ le-v){,-v') 

(a" + €)« "^ (fc" + £)» ■*■ (c'+e)» (e + o"; (e H- 6') (e + c") ' 

. (^\" , /^^^V , /<?6Y_ (e + a')(e + y)(6 + c') 

• • W "^Uyy "^ \dz) ~* (6-t/) (e-i;') ' 

* • W Uy W/ (e -«;)(«- w') * 

(t> — 1»; (v -e) (e — v) I ' da afy 

The equation V*F= is thus transformed into 
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or {v'-v) {(e + a«)(€ + 6»)(6 + c»)}i^1V 

+ (e-w') {(„ + o')(v+5')(v + c')}4~lV 

+ („_e) ^{(„'+„«)(v'+J«)(„'+c')}i^]V. 

5. A class of integrals of this equation, presenting a close 
analogy to spherical harmonic functions, may be investigated 
in the following manner. Suppose -E to be a function of e, 
satisfying the equation 

r{(€ + a')(€ + 5»)(€ + c«)}4^1V=(m€ + r)^, 

m and r being any constants. 

Then, if iff and H' be the forms which this function 
assumes when v and v are respectively substituted for e, 
the equation V'F= wiU be satisfied by V^EHH'. 

6. We will first investigate the form of the function 
denoted by E, on the supposition that ^ is a rational integral 
function of € of the degree n, represented by 

We see that 

+ + P,} 

= nr(n- 1) (c + a") (« + 6')(e + |6"-' + (n- 2)^,e-' 

(«-2)(n-3) ) 

+ 172 P*^ +-"+i>«-.| 



(n-1)(re-2') 



n-3 



F.H. 8 



..,}]. 
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Hence writing 

we see that 

= (me + r) je" + np/^ + "^."^^^ !>,£"- + ... +;),| . 

Hence, equating coefficienta of like powers of e, we get 
« f » + 2) = "»' 
n [(n - 1) {(n - 2)^, + S/J + 1 {(« - 1) ;). + 2/,}] = nmp, + r, 

||(^,-l)|-2)^^^2(«-l)/^.+/.}] 



^r 



w (n -- 1) 



or, as they may be more simply written, 
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n |(« - 1) ^i - 2) p* + Snf\ = nmp, + r, 

w (w — 1) 
= 1.2 ^ i^a + ^^i 



It thus appears that p^ is a rational function of r of the 
first degree, p^ of the second, p^ of the n**", and when the 
letters pp^g...^^ have been eliminated, the resulting equa- 
tion for the determination of r will be of the {n + 1)^ degree. 
Each of the letters p^, Pr^Pn ^^^ have one determinate 
value corresponding to each of these values of r; and we 

have seen that w = n ( ti + ^ j . There will therefore be (ti + 1) 

values of E, each of which is a rational integral expression 
of the 7jI^ degree, n being any positive integer. 

7. But there will also be values of E, of the n* degree, 
of the form 

(e+6')4(6+Oi|e->+(«-l)?,e"-'+^-^:^^^j,6-»+...+j,..}. 

We thus obtain. 

{(e + o«)(e + 5')(6+c')}*^ 

= (e + o»)i (e + 6')(e + <?){n - 1) je"-* + (» - 2) q^^ 



(n -2)(w-3) _^ 
1.2 



+ '" /'q "' ?.«■" + - + 2--.|; 

S— 2 
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+ (€ + a*) J (e + 6*)} (n - 1) |e"^ + (« - 2) q/^ 

(n-2)(«-3) 1 

+ 172"^^ — ?i« +-" + ?«-t| 

+ (e+o»)* (e+i*) (e+0 (n-1) (n-2) |€-»+(«-3) j,6-* 

Hence 
II {e + b^ (£ + (>•) +(€ + o') (e + c*) + (e + o')(6 + J*)| 

+ (e+a») (6+6*) (e + (n-l)(»i-2)|e^ + (n-3)j.6"-* 
= (me + r) je""' + (n - 1) jj€ 



+ ^ 7 w ' ?,€""+..•+ ff-x} ; 



(«-l)(«-2) _^ 
1.2 



.•.(«-l)(| + n-2) = m, 

(«-l){2a»+|(6» + cO + |(n-2)?,| 

+ (n-l)(n-2){o*+6» + c?' + (n-3)g^,} = (n-l)m2', + r, 
(n - 1) {(^+ aV + a'y) j,^ + (n - 2) a*JVj»^} =rj^.. 
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By a similar process to that applied above, we shall find 
that T is determined by an equation of the w*^ degree, and 

that m=(w — 1) f w — ^j, and that each of the letters j,, 

?«---?n-i ^s a rational function of r. Thus, there will be n 
solutions of the form 

There will also be n solutions of a similar form, in which 
the factors (e + c*)* (e + a*)*, (e + a*)* (e + 6*)* are respectively 
involved. Hence, the total number of solutions of the ri^ 
degree will be 4n + 1. 

8. We may now investigate the number of solutions of 

the degree n-^^, n being any positive integer. These will 

be of the following forms : three obtained by multipl3ring a 
rational integral function of e of the degree w by (e + a*)*, 
(e+6')*, (e + c*)*, respectively, and one by multiplying a 
rational integi-al function of e of the degree n — 1 by the 
product 

{(€ + a«) (€ + &•) (€ + 01*. 

An exactly similar process to that applied above will 
shew us that there will be n + 1 solutions of each of the 
first three kinds, and n of the fourth. Hence .the total number 
of such solutions will be S(n+l) + n, or 4m, + 3, that is 



H) 



+ 1. 



To sum up these results, we may say that the total 
number of solutions of the n^ degree is 4m + 1, n denoting 
either a positive integer, or a fraction with an odd numerator, 
and denominator 2. 

Similar forms being obtained for H, H', we may proceed 
to transform the expression EHH' into a function of a?, y, z. 

9. Consider first the case in which 
T, . «i-i ri(n— 1) -_- 
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Write this under the form 

^=(€-tt)i) (€-o)J...(6-a>J. 

Then £"=(1; — co^) (v — ©J ... (v— ©J, 

^' = (v'-a,0(./-o,.)...(v'-o,.). 

Hence 

^ir5'= (e-o),) (u-o),) (i;'-co,)...(€ -a)J (t;-a)„) (v - o>J. 

Now we have shewn (see Art. 4 of the present Chapter) 
that (e — ft»j) (i; — ©J (i;' — ©J 

Each of the factors of EHH' being similarly transformed, 
we see that EHH' is equal to the continued product of all 
expressions of the form 

(o> + a»)(a>+5')(a, + c')f^ + ^+-j-^^^ iV 

^ ^ ^ '^ ^ \ar-\-o) o+G) (T + cD /' 

the several values of © being the roots of the equation 

« . n-1 . ^ ("^ -■ 1) «-o 

As this equation has been already shewn to have (w 4- 1) 
distinct forms, we obtain {n + 1) distinct solutions of the 
equation V*F=0, each solution being the product of n 
expressions of the form 

That is, there will be w + 1 independent solutions of the 
degree 2n in x, y, z, each involving only even powers of the 
variables. 

10. To complete the investigation of the number of solu- 
tions of the degree 2n, let us next consider the case in which E 

=(6+y)i(e:l-c')i{e-'+(n-l);),e-^+ ^^\)^^-\ e""+...tp.-j. 
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The object here will be to transform the product 
(€ + 5^)* (i; + y)i (i;' + Vf (€ + O* {v + O* W + (^\ 

since the other factors will, as already shewn, give rise to the 
product of n — 1 expressions of the form 



a? , y , z 



,8 



+ T/r-+-rT— -1. 



a* + CO 6* + CO c* + CO 

Now, by comparison of the value of a? given in Art. 4, 
we see that 

Hence, we obtain a system of solutions of the form of 
the product of {n — 1) expressions of the form 



^+r.^ + 



a'+o) V+o> c^ + <o 



•1, 



multiplied by yz. Of these there will be w, and an equal 
number of solutions in which zx^ xy^ respectively, take the 
place of yz. 

Thus, there will be 4^1 + 1 solutions of the degree 2n in 
the variables of which n + 1 are each the product of n 
expressions of the form 

a« + co^6« + co^c* + co ' 
n are each the product of (n— 1) such expressions, multiplied 

n ... ... ... zx, 

n ... ... ... xy. 

11. We may next proceed to consider the solutions of the 
degree 2n + 1 in the variables a?, y, z. 

Consider first the case in which 
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Here the product {e + a*)* (v + a*)* (v + a*)* will, as just 
shewn, give rise to a factor x in the product EHH'. 

Hence we obtain a system of solutions each of which is 
the product of n expressions of the form 






a*+tt) i' + c» c^-fo) 



-1, 



multiplied by x. Of these there will be w + 1, and an equal 
number of solutions in which y, «, respectively take the 
place of the factor x. 

Lastly, in the case in which 

we see that in EHH* the product 

{eWf {v'^a'f (v'+a*)* (e+i*)* (iz+J*)* (u' + J^* (e+c')* 

(i;+c')*(v'+c')i 
will give rise to a factor xyz. 

Hence we obtain a system of solutions each of which is 
the product of [n — 1) expressions of the form 

^ . f ._l_ 1 



a* + a) 6* + a) c^ + o) 

multiplied by ccyz. Of these there will be n. 

Thus there will be 4n + 3 solutions of the degree 2n + I 
in the variables, of which 

(n + 1) are each the product of n expressions of the form 

-i-; 1- T^ — I- -Tj 1 multiplied by a?, 

(n + 1) are each the product of n such expressions, multiplied 
(n + 1) ... ... ... ... z^ 
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n are each the product of (n — 1) such expressions, multi- 
plied by xyz. 

12. Now an expression of the form G . EHH^, G being 
any arbitrary constant, is an admissible value of the potential 

/vt* «t* g>^ 

at any point within the shell -2+*^+ -5 = 1. But it is 

not admissible for the space without the shell, since it 
becomes infinite at an infinite distance. The factor which 
becomes infinite is clearly -B, and we have therefore to 
enquire whether any form, free from this objection, can be 
found for this factor. We shall find that forms exist, bearing 
the same relation to E that zonal harmonics of the second 
kind bear to those of the first 

Now considering the equation 

r{(6 + a»)(6 + y)(€ + 0}i^1V=(m€+r) R 

which we suppose to be satisfied by putting Z7= E, we see 
that, since it is of the second order, it must admit of another 

particular integral. To find this, substitute for U, E Ivde, 

we then have 

[{(6 + a')(e + J')(e+(0}i^jj7 

= [{(e+ a») (e + b') (e + c*)}* ^1 E.jvde 

+ {(e + a") (e + 6») (e + c^}i ^; 

= [{(6 + o*) (e + 6») (e + c»)}4 ff E . jvde 

+ 1 {(e + i*) (e + c*) + (e + c') (e + a») + (6 + o') (6 + J*)} £» 
+ (e+a«)(e + i')(e + 0(^.z, + £:|). 

j by Google 
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"Now, wnce hj Bupposition, the eqxraticm for the determi- 
nation of [7 is satisfied by putting [7= E, it follows that 

when JS jvcU is substituted for U, the terms involving ivde 

will cancel each other, and the equation for the determina- 
tion of V will be reduced to 

Idv 2dE If 1 1 ,1\ 

whence logv + 2 log E+log{{e + o*) (e + 6*) (e + c")}* 

* =logt;(,+ 21og^,+ loga6(J, 
», and E^ being the values of v and E, COTresponding to e = 0. 

Hence v = v^-^ — . ; 

'E'{(e + a'){e + V){€ + d'}]i' 

We may therefore take, as a value of the potential at 
any external point, 

V^v^E^ahcEHE'C ^ ^. 

For this obviously vanishes when e= oo . It remains so 
to determine v^ that this value shall, at the surface of the 
ellipsoid, be equal to the value G. EHH\ already assumed 
for an internal point. This gives 

Hence, putting v^ . E^ . dbc = F^, we see that to the value 
of the potential 

V EHH' r — ~ 

Jo Je^{(e + a*}(€ + i*J(6 + c')}4' 
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for any internal point, corresponds the value 



v^eeh: 



'^[ 



d€ 



for any external point. 

13. We proceed to investigate the law of distribution of 
density of attracting matter over the surface of the ellipsoid, 
corresponding to such a distribution of potential. 

Now, generally, if Bn be the thickness of a shell, p its 
volume density, the difference between the normal compo- 
nents of the attraction of the shell on two particles, situated 
close to the shell, on the same normal, one within and the 
other without will be im-pSn. This is the attraction of the 
shell on the outer particle, minus the attraction on the inner 
particle. 

But the normal component of the attraction on the outer 

particle estimated inwards is — -i— , 

And, if F' denote the potential of the shell on an in- 
ternal particle, the normal component of the attraction on 

dV 
it estimated inwards is — -,— . 

an 

Hence 47rpow = -^ -y— . 

'^ an an 

dV^dVdxdrdyd/Fda 
dn dx dn dy dn dz dn' 

doc 
And T- is the cosine of the inclination of the normal at 
dn 

the point a^y, z to the axis of x, and is therefore generally 
equal to e a , > ^ denoting the perpendicular from the 
centre on the tangent plane to the surface 

^ ^ f _^ ^ 



a'^^e V-^-e c^-^e 



= 1. 
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And we have shewn that 



(a''-6»)(a*-c*) ' 



t 2dx 1 

whence - -^ = -a^— . 

OS adx 
a* + 6 oe ' 

dx dx 

' ' dn" de' 

' ■ dn \dx de dy de dz de) de' 

Similarly ^' = 2e^. 

dn de 

Now r=V^.EHH'^ ^ 



o^'{(6 + a')(e + i»)(e,+c')l*' 

• ^^=F HH'^r ^ 

" de "' deJo E*{{e + a*) (e + i') (e + c*))*' 

therefore, generally, 

— = F HH'—r ^ 

de •* dej, £- {(a' + e) {b' + e) {c* + e)}* 

__ TT EHH' - — 

But, when the attracted particle is in the immediate 
neighbourhood of the surface, € = 0. Hence, the first line 
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dV . . . dV 

of the value of -r- becomes identical with the value of -r- » 
cLe de 

and we have 

dV dV_yEn' 1 

de de ' E, abc' 

E^ denoting the value which E assumes, when 6 = 0. 

Hence, . iirpBn = 2e T^ — ^ r- . 

But Bn, being the thickness of the shell, is proportional to 
e, and we may therefore write F" = ^ > ^ being the thick- 
ness of the shell at the extremity of the greatest axis ; 

r, a 1 HE' 
• ' f"^2wBaabc E, ' 

and this is .proportional to the value of F corresponding to 
any specified value of €, since HE' is the only variable 
factor in either. 

Hence functions of the kind which we are now considering 
possess a property analogous to that of Spherical Harmonics 
quoted at the beginning of this Chapter. On account of 
this property, we propose to call them Ellipsoidal Harmonics, 
and shall distinguish them, when necessary, into surface and 
solid harmonics, in the same manner as spherical harmonics 
are distinguished. They are commonly known as Lamp's 
Functions, having been fully discussed by him in his Lefons. 
The equivalent expressions in terms of x, y, z have been con- 
sidered by Green in his Memoir mentioned at the beginning 
of this chapter, and for this reason Professor Cayley in his 
"Memoir on Prepotentials," read before the Royal Society 
on June 10, 1875, calls them " Greenians." 

We may observe that the factor 

47r Ba ahc 

1 11 

is equal to . , ^ , and therefore also to z ^ or -. — r^ i 
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Hence^ it is equal to 



47r ' 

-^ {bcSa + ca£b + abBc) 



volume of shell ' 
and the potential at any internal point 

=s i volume of shell X EE^.o \ '- r» 

and the potential at any external point 

=: J volume of shell xZK.p T — 5-; 

^here for p must be substituted its value in terms of v and v. 

14. We will next prove that if F^, F, be two different 
ellipsoidal harmonics, dS an element of the surface of the 

ellipsoid, lie F^F, Jfi^=0, the integration being extended all 

over the surface. 

We have generally 

///(F.7T.-r,vr.)&ij,&=//(r,5-F.^)M 

And throughout the space comprised within the limits of 
integration, V* Fj = 0, V F^ = 0. Hence 

Now it has been shewn already that Fj, F, are each of 
the form EHECy where ^ is a function of e only> H the same 
function of u, H' of v. We may therefore write 

andsimilarly F. =/, (e)/. (v)y; (u'). 
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Hence r^=Vv£^ 

Hence •'i ^^ - ^ »'. y^ (^^ . 

Now, all over the surface, 6 = 0. Hence 

Hence, unless } \rl — i u\\ ^^> which cannot happen 

unless the functions denoted by f^ and f^ are identical*, or 
only differ by a numerical factor, we must have 



jjeV.V^dS^O. 



Now e is proportional to the thickness of the shell at 
any point. Calling this thickness 8e, we have therefore 



IS' 



SeV,V^d8=0. 



Hence, adding together the results obtained by integrating 
successively over a continuous series of such surfaces, we get 



jjjv^V^dxdydz==0; 



F,, V^ now representing solid ellipsoidal harmonics, and the 
integration extending throughout the whole space comprised 
within the ellipsoid. 

* This may be shewn more rigorously by integrating through the 
space bounded by two confocal ellipsoids, defined by the yalaes X and /*> of e. 
We then get, as in the text, 

Now the factor within { } cannot vanish for all valnes «f X and jt^ nnless the 
functions devoted by j\ and /, be identical, or only differ by a numerical 
factor. 
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15. It will be well to transform the expression 



//" 



eV^r.dS 



to its equivj^ent, in terms of v, v. 

For this purpose we observe that if ds, ds be elements of 
the two lines of curvature through any point of the ellipsoid, 
dS^dsds'. 

Now, 
d«* is the value of da? + dy^^dz* when € and v are constant, 
ds^ ... ... ... 6 and V 

therefore if e and v do not vary, 

2dx ^ dv 

»•. do? =0-5 d^* 

2 a* + 1; 

Similarly dy^l^Jv, dz=^\^^dv'. 



Again, differentiating with respect to a> the expression 

a? «/» «" . 
-u- h TT^ h -r 1, we get 



obtained for " — ^ ' 



a? f g' _ (v-ft)) (v -fi)) 
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therefore, putting w = v, 

«>* , y , -^ ^ (v-v)ie-v ) , 

A similar expression holding for ds* we get 

jei__iL (u'-t;)*(e-v)(f-v) , „ 

16 (a'+v) (6'+u) (c'+u) (a'+u') (6»+t;') (c'+w') •" " 

A • 1 _ 0^ , y* . ^* 

Again, ^-(a.^gy + (j. + g).+ (7qr^ 

~(a»+6)(6'' + 6)(c'+e)' 
writing c for a in the expression above ; 

It has been shewn that, integrating all over the surface, 
the limits of u are - c* and — i^ those of v , — i* and — a*. 

Hence, F^, F,, denoting two different ellipsoidal har- 
monics 

/-*• /-*• VJ,{v'-v)dvdv' .^ 

J -f ; -»« {(a'+v) (6'+w) (c'+w) (a»+ v') (6*+i/') (c'+u')]* " 

The value of the expression 1 1 1 V*dxdyda, or its equiva- 
lent 

^j, r^ r^' r(v--v)^i;e?i;- 

J-6^J-a^{(a»+i;)(6»+i;)(c«+v)(a^+i;')(&*+u')(c'+t;'j]*' 

in any particular case, is most conveniently obtained by 
expressing F as a function of a?, y, z. 

F. H. 9 
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16. Before proceeding further with the discussion of ellip- 
soidal harmonics in general, we will consider the special case 
in which the ellipsoid is one of revolution. We must enquire 
what modification this will introduce in the quantities which 
we have denoted by a, jS, 7, viz. 

^r d± 



Jv (a' + 






d^|r 



and in the differential equation 

We will first suppose the axis of revolution to be the 
greatest axis of the ellipsoid, which is equivalent to supposing 
6* = cl To transform a and 7, put a* + 'ilr = ^, a* -f e = rf, 
a* + u = 6)* ; we then obtain 

y Jo ^-a'+6^-(a«-i0^^^ («'-**)* + «' 
To transform ^, we must proceed as follows. 

Put i/r = - c'cos*tsr — ft^sin'-BT, v = — c* cos' <^ — 6' sin* <^, 
we then get generally 

6' + ^= (6»- c') cosV, c' + ^= (c*- 6") sinV ; 
dy^=^2 (c^ — 6*) cos -BT sin w d/^ ; 
. g_ 2 fo Jtsr ^2\g.0 
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Hence. ^ = -1 (,._„«+ j.) | . 

d _ (rt'-&')i d 

Also, e^tf — a*, t>' = ©* — «*, u = — 6',and our diflferential 
equation becomes 

(a«-6'-«»)|(,«_a» + 5')|}V 

-(a,'-i,')(a«-&')^=0, 
or («» - a* + 6») j(7,» - a» + J*) j-l' F 

_(a«-6')(,»-a,')^=0. 

This equation may be satisfied in the following ways. 

First, in a manner altogether independent of <f>, by sup- 
posing V to be the product of a function of rj and the same 
function of o), this function, which we will for the present 
denote by/ (17) or /(a)), being determined by the equation 



dri 
— 



{(«'-«' + 6*) ^}/H = 'n/W- 



9-2 
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Secondly, by supposing ^-, a constant m\iltiple o{ V, 
= — cr*F, suppose. 

Our equation may then be written 

- ff' (a* - V) {(«* -a* + V)- {rf - o» + h*)\ 7= 0, 

which may be satisfied by supposing the factor of F inde- 
pendent of ^ to be of the form F(;ii) F(ai), where 

](,»_ o» + 5») ^Fir,) -<^{a*- V) F{fi) = m (ij'-a'H- J») ^(i,). 

ko*-a* + V) ^|V<«) -<7» (a»-6') i?'(a))=«i(a>*-a»+J') !'(«»). 

The factor involving ^ will be of the form 

A cos o-^ +.6 sin a^. 
Now, returning to the equation 

||(,._aVJ')|}/(,) = V(,), 

we see that, supposing the index of the highest power of 17 
involved ^f(ji) to be i, we must have «i = i (« + 1). 

Now, it will be observed that ri may have any value 
however great, but that a>\ which is equal to c? + v\ must 
lie between a* — 6' and 0. Hence, putting ©^ = (a* — i') /a', 
where /a* must lie between and 1, we get 
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Hence this equation is satisfied by /{fa*— 6')*/i,} = (7P„ 
C being a constant ; and supposing C = 1 we obtain the 
following series of values for/ (o)), 

t = 0. /(a,) = l, 

* = 2' •^(")° 2(a'-y) ' 

Exactly similar expressions may be obtained for f(rf), and 
these, when the attraction of ellipsoids is considered, will 
apply to all points within the ellipsoid. But they will be 
inadmissible for external points, since rj is susceptible of in- 
definite increase. 

The form of integral to be adopted in this case will be 
obtained by taking the other solution of the differential 
equation for the determination of f{f))t i.e» the zonal har- 
monic of the second kind, which is of the form Q^ < — — —A , 
where 

Or, putting 17' = (a* — J*) i;*, ^ = (a* — h^ \*, we may write 

17. We may now consider what is the meaning of the 
quantities denoted by 17 and &). They are the values of Sr 
which satisfy the equation 



Digitized by VjOOQ IC 



134 ELLIPSOIDAL AND SPHEROIDAL HARMONICS. 

and are therefore the semi-axes of revolution of the surfaces 
confocal with the given ellipsoid, which pass through the 
point Xy y, z. One of these surfaces is an ellipsoid, and 
its semi-axis is 17. The other is an hyperboloid of two sheets 
whose semi-axis is a>. 

Now, if ^ be the eccentric angle of the point x, y, z, 
measured from the axis of revolution, we shall have 

cf^jf cos* 0. 

But also, since 17', w', are the two values of ^ which 
satisty the equation of the surface, 

Hence cd* = (a* — V) cos' 6, 

and we have already put 

whence the quantity which we have already denoted by fi 
is found to be the cosine of the eccentric angle of the point 
X, y, z considered with reference to the ellipsoid confocal 
with the given one, passing through the point a?, y, z. We 
have thus a method of completely representing the potential 
of an ellipsoid of revolution for any distribution of density 
symmetrical about its axis, by means ^of the axis of revo- 
lution of the confocal ellipsoid passing through the point 
at which the potential is required, and the eccentric angle 
of the point with reference to the confocal ellipsoid. For 
any such distribution can be expressed, precisely as in the 
case of a sphere, by a series of zonal harmonic functions of 
the eccentric angle. 

18. When the distribution is not symmetrical, we must 
have recourse to the fonn of solution which involves the factor 
A cos cr<^ + jB sin <r<^. It will be seen that, supposing F to 
represent a function of the degree i, and putting m^i (t+1), 
the equation which determines F{fo) is of exactly the same 
form as that for a tesseral spherical harmonia For -^(9;), if 
the point be within the ellipsoid, we adopt the same form, 
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if without, representing the tesseral spherical haxmonic by 

^Z*^) I — ^— il > or T^^ {v), we adopt the form 
l(a - Vy) 

19. It may be interesting to trace the connexion of sphe- 
rical harmonics with the functions just considered This may 
be effected by putting i' = a". We see then that rj will become 
equal to the radius of the concentric sphere passing through 
the point, and rj^ — a^4- V will become equal to if. Hence 
the equation for the determination of /(i;) will become 

which is satisfied by putting /{rj) =7j\ or 7f^**^\ The former 
solution is adapted to the case of an internal, the latter to 
that of an external point. 

With regard to /(©), it will be seen that the confocal 
hyperboloid becomes a cone, and therefore o) becomes inde- 
finitely small. But fi, which is equal to — ^ —.y remains 

finite, being in fact equal to - or cos 0. Hence /(/t*) becomes 

the zonal spherical harmonic. 

Again, the tesseral equations, for the determination of 
^(17), jF (ft)), become 

which are satisfied by F{r)) =1;* or if^^*^K 

And, wilting for (o\ (a'-i*) fi\ we have, putting F((o) =x(m)> 

which gives x W = ^i^'^Kf^)- 
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20. We will next consider the case in which the axis of 
revolution is the least axis of the ellipsoid, which is equi- 
valent to supposing a^ = b\ To transform a and fi, put 
c^ + ir = 0\ c* + € = iy*, c' + v = a)*, we thus obtain 

To transform y, we must proceed as follows : 

Put ■^ = — a*8in*«r — J'co8*«r, «' = — a'sin*^— i*cos*0, 
we then get, generally, 

a' + ir={a'-b*)co%*w, fc» + ^=-.(o*-fc*)8in»w, 

c*+'^=c*— o*sin*^— J*co8*^, rf^s— 2 (a*— fc*) sinw cos«r dar. 

Hence 

~ = 2 I r= ^^- — r if O = . 

^ i« (o» ainV + i' cos'isr - c")* (a» - c')* 

Hence, |. = - ^ (a«_c» + ,«)|. 



also, 



W = 0)* — C*, 

V = — a , 



and our differential equation becomes 

(a«-(?+a,»)|ca«-c* + ,')|}V 

- (a»-c'+iy') |(a»-c' + a,') ^|V 
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We will first consider how this equation may be satisfied 
by values of F independent of <f>. 

We may then suppose Fto be the product of a function 
of V, and the same function of o), this function, which we will 
suppose to be of the degree i, being determined by the 
equation ' 

^|(a«-c«+a,*)^|/(«,)=.-(t+i)/(a,). 

On comparing this with the ordinary differential eqtia- 
tion for a zonal harmonic, it will be seen that, on account 
of a' being greater than c*, the signs of the several terms in 
the series for /(rj) will be all the same, instead of being 
alternately positive and negative. We shall thus have 



{o»-c')4' 






and generally 



We will denote the general value of/ (17) hjPi \r^ — ^-rf , 
or, writing rj=:{a*^ c*)*^, by Pi{v). 
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For external points, we must adopt for /(17) a function 
which we will represent by qA t » A i » ^^ liW» '^tich will 
be equal to 

It is clear that /(co) may be expressed in exactly the 
same way. But it will be remembered that 97' and «* are 
the two values of y which satisfy the equation 

Hence 17, as before, is the semi-axis of revolution of the 
confocal ellipsoid passing through the point {x, y, z). But 
rfcci* = — (a* — c*) «', an essentially negative quantity, since 
a' is greater than c*. Hence co' is essentially negative. » Now, 
if ^ be the eccentric angle of the point {x, t/, z) measured 
from the axis of revolution, we have «' = if cos'd. Hence 

,7V = - (a* - c^ 17* cos*^, 

and therefore co* = — (a* — (?) cos' 

= — (a' - c") /A*, suppose. 

Hence the equation for the determination of /(©) assumes 
the form 

the ordinary equation for a zonal spherical harmonic. Hence 
we may write 

fi being the cosine of the eccentric angle of the point a?, y, z, 
considered with reference to the confocal ellipsoid passing 
through it. 
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21. We have thus discussed the form of the potential, 
corresponding to a distribution of attracting matter, sym- 
metrical about the axis. When the distribution is not 
symmetrical, but involves ^ in the form A cos <r<f> + B sin <7</>, 
we replace, as before, P, (/a) by T/*'^ (/a), and Pi(fi) by a 
function t^^'^^^v) determined by the equation 

<.»(.) = (l + ;<')^^ ft M, 
and 2. (V) by ,(.)(.) /J ^^^|^-_. 

22. As an application of these formulae, consider the fol- 
lowing question. 

Attracting matter is distributed over the shell whose 

surface is represented by the equation -^+ j^ =1, so 

that its volume density at any point is P, (fi), fi being the 
cosine of the eccentric angle, measured from the axis of 
revolution ; required to determine the potential at any 
point, external or internal. 

The potential at any internal point will be of the form 

CF,MP,(v) (1), 

and at an external point, of the form 

(rp,MQ,{v) (2), 

where (a* — 5*)* i/ = the semi-axis of the figure of the con- 
focal ellipsoid of revolution passing through the point (fi, v). 

Now the expressions (1) and (2) must be equal at the 
surface of the ellipsoid, where v = — — — . . 

Hence 
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But generally 



Hence 



We may therefore, putting (7' = u4P, J ^ t^m ' ^^ 

and we thus express the potentials as follows : 

AP^ Qi) P^ (y) Q^ \-— — — -, h at an internal point, 
({a —o)V 

AP^ (ji) Q^ (v) P, ] — ^j- at an external point. 

Or, substituting for Q^ its value in terms of P,, 

at an internal point, 

at an external point. 

Now, to determine A, we have, Sa being the thickness 
of the shell at the extremity of the axis of revolution, 
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'^ 4iTrSa.i]\drj dr) Jy,-a 



4iTr'Sad'-b*\dv "" dv r^'-f^ 



- i^ kir^* ^^' ^>^' 1(^^41 ^' i(^^4j 



la -4 _ . . 1 



a" -6' 



.-1 



= i:;;ws:^^«0»)- 



4ir 6'3a' 
Hence, if p = P, (/*), we obtain 

a 
And we thus obtain 



K = 47r 6S6 P, 0*) P (v) P, ( " . 1 r ^ 

= 4^5S6P.(m)P.W«.{^j^:^4}. 

r.=4.^P,(.)P.{^-^,}p,WJ;=^^^ 
= 4.6SJP,W(2.(.)P,|^^}. 



If the shell be represented by the equation 
it may be shewn in a similar manner that we shall have 
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r, = 4i7radaP, (fi) q^{p) p, L,^^,J • 

23. We may apply this result to the discussion of the 
following ptoblem. 

If the potential of a shell in the form of an eUipsoid of 
revolution about the greatest point be inversely proportional 
to the distance from one focus, find the potential at any 
internal point, and the density. 

If the potential at P be inversely proportional to the 
distance from one focus 8, and H be the other focus, we have, 

HP +8P-= 27,, HP^8P^%o, 

/. SP = v-^' 

Hence if M be the mass of the shell, V^ the potential at 
any external point,. 

' tf — m 

M 1 



(o'-6»)ii'-/* 

Now, by what has just been seen, the internal potential, 
corresponding to P^ (ji) Q^ [v), is 

^'|(a«-6')4|. 
Hence, if V^ be the potential at any internal point. 

Digitized by VjOOQ IC 



ELLIPSOIDAL AND SPHEBOIDAL HABMOKICS. . 143 

^' = (^4 ^ (2i + 1) -^l^^l P. 0*) P, W. 
And the volume density corresponding to P, (ji) Q^ {v) is 



^mp. 



.[ ^ r 



Hence the density corresponding to the present distri- 
bution is 

If F, had varied inversely as HP, we should have had 

and our results would have been obtained from the foregoing 
by changing the sign of q), and therefore of fi. 

24. Now, by adding these results together, we obtain 
the distributions of density, and internal potential, corre- 
sponding to 

or, in geometrical language,. 

= M multiplied by the axis of revolution of the confocal 
ellipsoid, and divided by the square on the conjugate semi- 
diameter. We may express this by saying that the potential 
at any point on the ellipsoid is inversely proportional to the 



Digitized by VjOOQ IC 



144 ELLIPSOIDAL AND SPHEBOIDAL HABMONIC& 

square on the conjugate semi-diameter, or directly as the 
square on the perpendicular on the tangent plane. 

Corresponding to this, we shall have, writing 2k for t, 
since only even values of i will be retained, 



(a*-i')*" 



^A{a*lb')i 









k being 0, or any positive integer. 

Again, subtracting these results we get 
^ if Jf ,, 2tt> 

= M multiplied by the distance from the equatoreal plane, 
and divided by the square on the conjugate semi-diameter. 

This gives, writing 2k + 1 for i, 









25. In attempting to discuss the problem analogous to 
this for an ellipsoid of revolution about its least axis, we see 
that since its foci are imaginary, the first problem would re- 
present no real distribution. But if we suppose the external 
potential to be the sum or difference of two expressions, each 
inversely proportional to the distance from one focus, we 
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obtain a real distribntion of potential — in tbe ftrsi case 
inversely proportional i to the square on the conjugate 
semi-diameter, in the latter Varying as the quotient of the 
distance from the equatoreal plane by the square on -the 
conjugate semi-diameter. 

It wUl be found, by a process exactly similar to that just 
adopted, that the distribution3 of intern^ potential, and 
density, respectively corresponding to these will be : 



In the first case 



g"V-c4 



F, = ^iS (4*4-1) 4±l£4p,C)^„(.). 



^"wh^l 






Je being 0, or any positive i^teger. 
In the second case 



^^= rJ^iS(4A;+8) ^'^y -/''i p^Mp^M 



f g ] 

'ICa'-c*)^. 






Pttfi- 



(a«-c')iJ 



Z; being 0, or any positive integer, 

26. We may now resume the consideration of the ellip- 
soid with three unequal axes, and may shew how, when the 
potential at every point of the surface of an ellipsoidal shell 
IS known, the functions which we are considering may be 
employed to determine its value at any internal or external 
point. We will begin by considering some special cases, 

F. H, 10 
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by Avhicli the general principles of tlie method may be made 
more intelligible, 

27f . First, suppose that the potential at every point of 
the surface of the ellipsoid is proportional to a? :== — ^ suppose, 

Ia this case, since x when substituted for V, satisfies the 

equation v*^= 0, we see that V^- will also be the potential 

at any internal point. But this value will not be admissible 
at external points, since x becomes infinite at an infinite 
distance. 



»-{^ 



Now, transforming to elliptic co-ordinates 
f (e+ffl')(v+a')(i;' + a«) U 

And the expression 
Vj, ( (6+a*)(i>fa')(v'+a') U f d^ . 



Jq I 



satisfies, as has already be€fn seen, the equation y*^=sO, is 

equal to Fo- at the surface of the ellipsoid, and vanishes 

at an infinite distance. This is therefore "^the value of the 
potential at any external point. It may of course be written 



< 



{f'\-a'){(Sr + a') if + b') (i^ + c")}** 



28. Next, suppose that the potential at every point of 
the surface is proportional to y«=Fo?-, suppose. In this 
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ca^, as in the last, we see that, since yz when substituted 
for V, satisfies the equation ^W^ 0, tne potential at any 

iijtemal point will be F. ^; while, substituting for y, z their 

values in terms of elliptic co-cnrdinates we obtain for the 
potential at any external point 

*« /• (t + 6') (t + f (t + a*) {f + &•) (t + <^)f 



■r, 



29. We will next consider the case in which the po- 

tential, at every point of the surface, varies as o?=^V^-i 

Suppose. This case materially differs from the two just con- 
sidered, for since a? does not, when substituted for F, satisfy 
the equation v*F= 0, the potential at internal points capnot 
in general be proportional to a?. We have therefore first to 
investigate a function of a?, y, «, or of €, v, v which shall 
satisfy the equation v"F=^0, shall not become infinite within 
the surface of the ellipsoid, and shall be equal to as* on its 
surface. 

Now we know that, generally 
(y + «d) (c" + 0)) a?+ (c*+ ©) (a» + o)) y" + (a* + ©) (6* + ©) «* 

-(a' + ii))(i* + €«)) (c' + «D) = (€-ii))(i;-G)) (i;'-6)). 

And, if '^^, 0^ bQ the two values of <o which satisfy the 
equation 

we see that 

and ^(e-e^{v-e;){v'-e^^Q. 

And, by properly determining the coefficients A^, A^, A^, 
it is possible to make 

A,+A,ie-e;)(y-exv'-0d+A»{e'e;){v-0Xv'-0:).A2) 

^ ^ whpn 6'cV + cW + aV^ - a»6V = 0. 
or ' " 

10—2 
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Hence, the expression (2) when A^, -4^, A^ are properly 
determined, will satisfy all the necessary conditions for an 
internal potential, and will therefore be the potential for 
every internal point 

Now, we hJtve in general 

-{a' + 0Xb' + O^)i(^ + O^^{e-0,)(y''ff^W-ff^) 
(b' + 0;){c\+0,)a? + ((f + 0;)(a' + 0;)f+{a'-h0,)(b' + 0;}s? 

and, over the surface 

6'c V + c'oy + 0*6 V - o*&V = 0. 

Hence, Sf being any quantity whatever, we have, all over 
the surface, 

and therefore, putting ^ = — a'. 

Hence, the right-hand member of this equation possesses 
all the necessary properties of an internal potential. It 
satisfies the general dififerential equation of the second order, 
'does not become infinite within the shell, and is proportional 
to x^ all over the surface. 

.We observe, by equation (1), that 
(&«+©)(c^+a)) + (c'+©) (a»+a>) + (a'+o) (&'+a))= 3(^,-0)) {0,-w) 
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•identically f and therefore, writing —a* for ai, 

(a»-5^(a«-c*)=.3(a*+.^^(a» + ^^. 

Hence, over tlie surface of the shell, 

S 2 

and we therefore have, for the internal potential, 

This is not admissible for external points, as it becomes 
i^finite at an in^nite distance. We must therefore substi- 
tute for the factor € — ^^ 

r, g>r d^ 

with a similar substitutioa for e — ^,, thus giving, for the 
external potential, 

; ^r §±_ 

' Jo {f-e,y {(t + a') (^ + 6') if + c')j^ 

^r df' 
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The distribution of density over the surface, correspoDd- 
ing to this distribution of potential, may be investigated by 
means of the formula 

_ 1 a (dV^ dV\ 
^ 2.irda\d€ dej.^^' 

o>r its equivalent in Art 13 of this Chapter. We thus find that 

■■ 1 g n r („-g,)(v'-g.) 



■r< 



cZ^ 



{v-e:)W^0,) . r d^ 

30. The investigation just given, of the potential at an 
external point of a distribution of matter giving rise to a 
potential proportional to a? all over the surface, has an in- 
teresting practical application* For the Earth may be re- 
garded as an ellipsoid of equilibrium (not necessarily with 
two of its axes equal) under the action of the mutual gravi- 
tation of its parts and of the centrifugal force. If, then, 
F denote the potential of the Earth at any point on or with- 
out its surface, and li the angular velocity of the Earth's 
rotation, we have, as the equation of its surface, regarded as 
a surface of equal pressure, 

••. F+ ^ ft* (a? + ^) =* a constant, 11 suppose. 

Hence, if a, 5, c denote the semi-axes of the Earth, we 
"havo, for the determination of F, the following conditions: 
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7*= at an iofinite distance, < (2), 

F= n -- 1 Il» (a? + y») when 

The term 11 will, as we know^ give rise to an external 
potential represented by 

nf ^ ^^ - T ^^ 

;. {(^+a«)(^+6')(f +c»)}i • Jo {(^+a«)(t+6')(^+c*)}i* 

1 1 ' 

The two terms — « ^*^> "" o ^V* "'^ P^® ™® *^ terms 

which may be deduced from the value of F, just given by 

' 1 1 . 

successively writing for Y^ -^-^ ^<if% ttnd . — ^ ft'6^ and (in 

the latter case) putting V for a' throughout^ We thus get , 

^°(°-i"''"'^^''}f i(^^v)(A)(»:fc-)ii : 

*i.{(t+o')(^+6')(t+0}* 6 U+^» J'+V . 
f>i («!- ^,) i. (Vr - ^,)'' ((t + a«)(t + &•) (V' + c^l* 

^ r rf^ n^ f a' y \ 

d^fr 



r ^^ 
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31, Any rational integral function F of x, y, z, which 
satisfies the equation y*F*= 0, can be expressed in a series 
of Ellipsoidal Harmonics of the degrees 0, 1, 2...t in a?, y, z. 
For if V be of the degree i, the number of terms in V will 
^^ (♦+l)(t + 2)(t + 3) ^ Now the condition v«F=0 is 

equivalent to the condition that a certain function of x, y, z 
of the de^ee t— 2, vanish^ identically, and this gives rise 

to ^ ^-^ conditions. Hence the number of inde- 



pendent constants in F is 

(f + l)(tH-2)(t+3) (f^l)f(t + l) 
6 6 ' 

or {i+ 1)\ And iJie number of ellipsoidal harmonics of the 

degrees 0, 1, 2...t in tc, y, z or of the degrees 0, ^, 1, ^...3 

in 6, v^ V, is, as shewn in Arts. 6 to 10 of this Chapter, 

• l + 8 + 5+...+2t + l, 

or (i + 1)*. Hence all the necessary conditions can be satis- 
fied. 

32. Again, suppose that attracting matter is distributed 
over the surface of an ellipsoidal sheU according to a law of 
density expressed by any rational integral function of the 
co-ordinates. Let the dimensions of the highest term in this 
expression be i, then by multiplying evenr term, except those 
of the dimensions % and » — 1 by a suitable power of 






we shall express the density by the sum of two rational inte- 
gral functions of a?, y, z of the degrees i, t — 1, respectively. 
The nmnber of terms in these will be 
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And any ellipsoidal surface harmonic of the degree ^ i — 2... 
^ ^1 y; ^> ^3^7> hj suitably intitnlucihg the £Eu;tor 

a? v" s? 

be expressed as a homogeneous function of rv, j/, z of the 
degree i ; also any such harmonics of the degree t— 1, t— 3... 
in X, y, z may be similarly expressed as a homogeneous 
function of x, y, z of the degree % — 1. And the total number 
of these expressions will, as just shewn, be (t + 1)*, hence by 
assigm'ng to them suitable coeflScients, any distribution of 
density according to a rational integral function of x, y, z 
may be expressed by a series of surface ellipsoidal harmonics, 
and the potential at any internal or external point by the 
corresponding series of solid ellipsoidal harmonics. 

33. Since any function of the co-ordinates of a point on 
the surface of a sphere may be expressed by means of a series 
of surface spherical harmonics, we may anticipate that any 
function of the elliptic co-ordinates v, v may be expressed by 
a series of surface ellipsoidal harmonics. No general proof, 
however, appears yet to have been given of this proposition. 
But, assuming such a development to be possible at all, it 
may be shewn, by the aid of the proposition proved in 
Art 15 of this Chapter, that it is possible in only one way, 
in exactly the same way as the corresponding proposition 
for a spherical surface is proved in Chcip. IV. Art. 11. 

The development may then be effected as follows. De- 
noting the several surfetce harmonics of the degree iiaXf y, z, 

or ^ in V, V, by the symbols F/i>, F/«>,-..F/2»-»'i), and by 

F{v, v) the expression to be developed, assume 

^-c;a>F,a)+... + c;wF/<^>+... 

Then multiplying by eF/'> and integrating all over the 
surface, we have 

jeF (v, v') F/'> dS = a/'> J^ (F/'>)* dS. 
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The values o{jeF{v, v')F/'> dS, and oije ( V^^ydB must 

be ascertained by introducing the rectangular co-ordinates 
a?, y, Zy or in any other "w».y which may be suitable for the 
particular case. The eoefficients denoted by C are thus 
determiifeed, tend fhe development effected^ 
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EXAMPLES. 



1. Prove that i^Oy^^P,-^P, + ^P,. ■ 

Why <»niiot (anff)* be expanded in a finite series of spherical 
harmonics} 

111 ^+«^« 

2. Prove that \ + ^P,+ ^ P, + j P',+ ... =log-— -^4 • 

3. Establish the equations 



(/^•-i)-i"=«/^.-«^.- 



diL 



i> 



nP. = (2«-l)/J'„-(»-l)P.... 

4. If fi = cos fl, prove that 

P.O.) = l-i(»H-l)sin'| + ... + (-ir^l^^(sin'|)%... 
and also that 

^*w=(-iy+(-ir*(*+i)<^o«"|+- 

(| w) l^-w \ 2/ 

5. Prove that, if a be greater tHan c, and * any odd 
integer greater than w>, 

6. Prove that /"V^*yrf/t=t(i+l> 
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7. Prove that when tt=* 1, -y-;;*= r: \v,/ . 



8. Prove that 



s 



IP P 

P P P 

••• ••• ••• ••• 

^i9 •*Hl ••• •* a«-l 



is a nnnierical multiple of 



9. Prove the following equation, giving any Laplace^s co- 
efficient in terms of the preceding one : 



P^ = pP,+njlP.dp + C, 



i 



where Op- [ifi + Ji - fi* Jl - fif' cos(w-ta') and C is zero if w be 
even, and 

(-l)'" 2.^.| ,^(^^l -)p .if«beodd. 

10. If i,jf i 1)6 three positive integers whose sum is even, 
prove, that 

" 1.3...(i+iSr-i~l) 1.3...(jg; + i~i~l) l.Z ...{i^j^k-l) 

1.3 ... (» + y + A?-l) i+J4-A?+l ' 
Hence deduce the expansion of F^Fj in a series of zonal 
har^I6nics. 

11. Express afy + y* + yz +y + z as a sum of spherical 
harmonics. - 

12. Find all the independent symmetrical complete harmpnics 
of the third degree and of the fifth negative degree. 

13.* .Matter is distributed in an indefinitely thin stratum over 
the surface of a sphere whose radius is unity, in such a manner 
that the quantity of matter laid on an element .{BS) of the surface 
is 8^(l+aaj+ Jy + c«+/aj* + 5ry' + A«*), 
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where x, y, st are rectangular co-ordinates of iie element 8*? re- 
ferred to the centre as origin, and a, h, c, /, g, h are constante. 
Find the value of the potential at any point, whether internal 
or external. 

14. If the radius of a sphere be r, and its law of density be 
p = ax + hy + cZf where the origin is at the centre, prove that its 

potential at an external point (f, ly, {) is t^Tm (^+ i^ + c£) vbere 

E is the distance of {(, rj, Q from the origin. 

15. Let a spherical portion of on infinite quiescent liquid be 
separated from the ^quid round it hj an infinitely thin flexible 
membrane^ and let this membrane be suddenly set in motion, 
every part of it in the direction of the radius and with velodty 
equsd to St, a harmonic function of position on the surface. Find 
the velocity produced at aoy external or internal point of the 
liquid. State the corresponding proposition in the theory of 
Attraction. 

16. Two circular rings of fine wire, whose masses are Jf and 
M, and radii a and a\ are placed with their centres at distances 
h, b\ from the origin. The lines joining the origin with the 
centimes are perpendicular to the planes of the tings, and are in- 
clined to one another at an angle 0, Shew that the potential of 
the one ring on the other is 



^^^^'{^^.KQn)/ 



and j^^ and Q^ are the same functions of V and a' and of cos d 
an d sin 6 respectively, and c is the greater of the two quantities 
N/a« + ^and V^+F*. 

17. A uniform circular wire, of radius a, charged with 
electricity of line-density e, surrounds an uninsulated concentric 
spherical conductor of radius c; prove that the electrical density 
at any point of the surface of the conductor is 

the pole of the plane of the wire being the pole of the harmonics. 
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18. Of two spherical conductors, one entirely surrounds the 
other. The iliner has a given potential, the outer is at the 
potential zero. The distance between their centres being so 
small 'that its square may be neglected, shew how to find the 
potential fit any point between th^ spheres^ 

19. If the equation of the bounding surface of a homo- 
geneous spheroid of ellipticity c be of the form 



r=«(i-|cP.), 



prove that the potential iat any external point wiU be 

where O and A are the equatoreal and polar moments of inertia 
of the body^ 

Hence prove that Fwill have the same value if the spheroid 
be heterogeneous, the surfaces of equal density difering fron^ 
spheres by jv harmonic of the second order, 

20. The equation i? = a(l+.ay) is that of the bounding 
surface of a homogeneous body, density unity, differing slightly 
in form and magnitude from a sphere of radius a; a is a 
small quantity the powers of which above the second may be 
neglected; and y ia a function of two co-ordinate angles, such 
that 

where T^ T .,.Z^, Z^,,. are Laplace's 'functions. Prove that 
the potential of the body's attraction on an external particle, 
the distance of which from the origin of co-ordinates is r, is 
given by the equation 

4iraVf o . w + 2 g'y ^ ) 

r(. 2r' 4w + 2r»* J 
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EXASiPLES, I5d 

21. If if be the mass of a tuiiform hemisplierical shell of 
radius c, prove that its potei^tial, at any point distaut r from 
the centre, willhQ 

It IMfl JLp'^ 

according as r is less or greater than c; the vertex of the hemi* 
sphere being at the point at which /x = 1* 

22. A solid is bounded by the plane of a^, and extends to 
infinity in all directions on the positive side of that plane. 
Every point within the circle a^ + y* = a*, z = is maintained at 
the untform temperature unity, and every point of the plane xt/ 
•without this circle at the uniform temperature 0. Prove that, 
when the temperature of the ^olid has become permanent, its 
value at a point distant r from the origin, and the line joining 
which to the origin is inclined at an angle d to the asps of z will 
be 

r 1 r» 1.3 r* 
P -P L 4.--P ^ P L + 

1.3...(2t^l) f^^ 
^ ^" ^^ 2.4...2i ^^^ a«+^ ■*■ ••* 
if r< a, and 

1 «» 1.3 a* ,,, 1.3...(2t-l) ^ a« 

2 ^ r* "274 ^»7* *** — " ^^ ^/ 2.4...2t -*+* r^* "* 

ifr>a. 

23. Prove that the potential of a circular ring of radius (?, 
whose density at any point is cos m^, c^ being the distance of the 
point measured along the ring j&om SQine fixed point, i^ 
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1 d;]^,^ 

(2.4.6...2WI dfiT f^* 



2ircoflm^(8ing)"| ^ ^ ^ o^ -:/::;;^^! 



1 <?"-?U <^. 

J '**2,4.6...(2w + 2) c&t" r-**"*" 

f '' . 1.3.5.,.(2A:^1) cTP,,^ <f^«* 

where r is greater tlian c. If r be iess than c, r and c mtust be 
y interchangecL 

24. A solid is bounded by two confocal ellipsoidal surfaces, and 
its ^Gi^^^y at any point P varies as the square on the perpendicular 
from the centre on the tangent plane to the confocal ellipsoid 
passing through P. Prove that the resultant attraction of such 
Ja solid on any point external to it or forming a piart of its mass 
is in the direction of the normal to the confocal ellipsoid passing 
through that point, and that the solid exercises no attraction on a 
point within its inner surface. 
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